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Abstract 

The resonant substructures of —)■ decays are studied with the Dalitz plot 

technique. In this study a data sample corresponding to an integrated luminosity 
of 3.0 fb“^ of pp collisions collected by the LHCb detector is used. The branching 
fraction of the —)• decay in the region m{D^7r^) >2.1 GeV/c^ is measured 

to be (8.46 ± 0.14 ± 0.29 ± 0.40) x 10“^, where the first uncertainty is statistical, 
the second is systematic and the last arises from the normalisation channel B^ —)• 
Zl*(2010)“7r'’'. The 7r'*“7r“ S-wave components are modelled with the Isobar and 
K-matrix formalisms. Results of the Dalitz plot analyses using both models are 
presented. A resonant structure at m(D^7T~) ~ 2.8GeV/c^ is confirmed and its 
spin-parity is determined for the first time as = 3“. The branching fraction, mass 
and width of this structure are determined together with those of the Z1 q( 2400)“ 
and 712(2460)“ resonances. The branching fractions of other B^ —)• decay 

components with —)• 7r+7r“ are also reported. Many of these branching fraction 
measurements are the most precise to date. The first observation of the decays 
B° DO/o(500), B^ DO/o(980), B^ D^p{U50), B^ 71|(2760)-7r+ and the 
first evidence of B^ —)• 7l‘^/o(2020) are presented. 
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1 Introduction 


The study of the Cabibbo-Kobayashi-Maskawa (CKM) mechanism 01 is a central 
topic in flavour physics. Accurate measurements of the various CKM matrix parameters 
through different processes provide sensitivity to new physics effects, by testing the global 
consistency of the Standard Model. Among them, the CKM angle is expressed in terms 
of the CKM matrix elements as eiTg{—VcdV*fj/VtdV^l). The most precise measurements have 
been obtained with the —)■ (cc)K^*^^ decays by BaBar [^, Belle and more recently 
by LHCb j^. The decay p] through the b —>■ cud transition has sensitivity 

to the CKM angle [3 (6}jl^ and to new physics effects 11-14 . 


The Dalitz plot analysis 15 of B^ —decays, with the iC+vr mode. 


is presented as the hrst step towards an alternative method to measnre the CKM angle 
/3. Two sets of results are given, where the 7r’''7r“ S-wave components are modelled 
with the Isobar 


50 ^ 


16 18 and K-matrix 19 formalisms 


TT 


have already been performed by Belle 


Dalitz p lot analyses of the decay 
and BaBar 


20 


21 


22 


Similar 


stndies for the charged B decays B~ —)■ have been pnblished by the B- 

factories 23 2^. The LHCb dataset offers a larger and almost pure signal sample. Feynman 
diagrams of the dominant tree level amplitudes contributing to the decay B^ — D^n~^7r~ 
are shown in Fig. [T} 



Figure 1: Examples of tree diagrams via b —)• cud transition to produce (a) resonances, (b) 

nonresonant three-body decay and (c) D^'k~ resonances. 


In addition to the interest for the CKM parameter measurements, the analysis of the 
Dalitz plot of the B^ —)■ decay is motivated by its rich resonant structure. The 

decay B^ —)■ contains information about excited D mesons decaying to Dvr, with 

natnral spin and parity = O’*", 1“, 2~^, ... A complementary Dalitz plot analysis of 
the decay B^ —>■ was recently pnblished by LHCb 25,26 , and constrains the 


^The inclusion of charge conjugate states is implied throughout the paper. 
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phenomenology of the D^K~ {D~j) and K~tt^ states. The spectrum of excited D mesons 
is predicted by theory 27,^ and contains the known states Z7*(2010), Z7 q( 2400), 172(2460), 
as well as other unknown states not yet fully explored. An extensive discussion on theory 
predictions for the cu, cd and cs mass spectra is provided in Refs. (^,29 . More recent 
measurements performed in inclusive decays by BaBar 


[30| and LHCb [^, have led to 
the observation of several new states: (2650), (2760), and Z7*(3000). However, their 

spin and parity are difficult to determine from inclusive studies. Orbitally excited D 
mesons have also been studied in semi-leptonic B decays (see a review in Ref. [^) with 
limited precision. These are of prime interest both in the extraction of the CKM parameter 
\Vcb\, where longstanding differences remain between exclusive and inclusive methods (see 
review in Ref. 32 ), and in recent studies of R —)■ D^*^TVr 
theoretical discussion (see, e.g.. Refs. 34 35] ). 


33 


which have generated much 


A measurement of the branching fraction of the decay R° —)■ is also presented. This 

study helps in understanding the effects of colour-suppression in B decays, which is due to 
the requirement that the colour quantum numbers of the quarks produced from the virtual 

Moreover, 


W boson must match those of the spectator quark to form a meson 36 -40 


using isospin symmetry to relate the decay amplitudes of R° —)■ D^p^, R° —)■ D~p^ 
and R+ —)■ effects of hnal state interactions (FSI) can be studied in those decays 

(see a review in Refs. 37,^). The previous measureme^ for the branching fraction 
of R° —)■ Z7°p° has limited precision, (3.2 ± 0.5) x 10“^ 
theoretical predictions that range from 1.7 to 3.4 x 10“^ 


21 


38 


and is in agreement with 


42 


Finally, a study of the tt+tt system is performed on a broad phase-space range in 
R° —)■ from 280MeV/c^ (~ 2m.n-) to 3.4GeV/c^ (~ m^o — m^o), which is much 

larger than that accessible in charmed meson decays such as —)■ 

B decays such as —)■ 


46 


49 


43-45 


The nature of the light scalar 


TT 


or in 
states 


below 1 GeV/c^ = O’*''*'), and in particular the /o(500) and /o(980) states, has been a 
longstanding debate (see, e.g.. Refs. |50| - |^ ). Popular interpretations include tetraquarks, 
meson-meson bound states (molecules), or some other mixtures, where the iso-singlets 
/o(500) and /o(980) can mix, therefore leading to a non-trivial nature (e.g. pure ss state) of 
the /o(980) and complicating the determination of the CKM phase (pg from —)■ 

decays |^[5^|^. In the tetraquark picture, the mixing angle, cUmix, between the /o(980) 

20° g 


and /o(500) states is predicted to be Iwmix 
average of the mass of the k meson 682 ± 29 MeV/c^ 32 


56] (recomputed with the latest 
). Other theory models based on 
QCD factorisation and its extensions ]57]|^ predict that the /o(500) and /o(980) mixing 

< 


angle (p-mix for the qq model is 20' 
of —>■ Jlp)'K^T\~ decays 

R° (R°) decay: (p^ix < 17° 
mixing in the B^ —)■ 


47 

49 

< 7 

.7° 


<Pmix ^ 45°. The LHCb experiment, in the study 
, has already set stringent upper bounds on pmix in 
at 90% CL. For the first time, the /o(500) — /o(980) 


TT 


decay, both in qq and tetraquark pictures, is studied. 

0 ., 


The analysis of the decay R° —?■ presented in this paper is based on a data 


sample corresponding to an integrated luminosity of 3.0 fb~^ of pp collision data collected 
with the LHCb detector. Approximately one third of the data was obtained during 2011 
when the collision centre-of-mass energy was ^/s = 7TeV and the rest during 2012 with 
/s = 8TeV. 
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The paper is organised as follows. A brief description of the LHCb detector as well 
as the reconstruction and simulation software is given in Sec. The selection of signal 
candidates and the ht to the candidate invariant mass distribution used to separate 
and to measure signal and background yields are described in Sec. An overview of the 
Dalitz plot analysis formalism is gAen in Sec. Details and results of the amplitude 
analysis hts are presented in Sec. In Sec. the measurement of the 
branching fraction is documented. The evaluation of systematic uncertainties is described 
in Sec. The results are given in Sec. and a summary concludes the paper in Sec. 


2 The LHCb detector 


The LHCb detector 

range 2 < r/ < 5, designed for the study of particles containing 6 or c quarks 


59 is a single-arm forward spectrometer covering the pseudorapidity 

The 


detector includes a high-precision tracking system consisting of a silicon-strip vertex 
detector surrounding the pp interaction region j^, a large-area silicon-strip detector 
located upstream of a dipole magnet with a bending power of about 4Tm, and three 
stations of silicon-strip detectors and straw drift tubes 61 placed downstream of the 


magnet. The tracking system provides a measurement of momentum, p, with a relative 
uncertainty that varies from 0.4% at low momentum to 0.6% at lOOGeV/c. The minimum 
distance of a track to a primary vertex, the impact parameter (IP), is measured with a 
resolution of (15 -|- 29 /pt) where px is the component of p transverse to the beam, 
in GeV/c. Different types of charged hadrons are distinguished using information from 
two ring-imaging Gherenkov detectors 62 . Photon, electron and hadron candidates are 


identihed by a calorimeter system consisting of scintillating-pad and preshower detectors, 
an electromagnetic calorimeter and a hadronic calorimeter. Muons are identihed by a 
system composed of alternating layers of iron and multiwire proportional chambers 


63 


The online event selection is performed by a trigger which consists of a hardware stage, 
based on information from the calorimeter and muon systems, followed by a software 
stage, which applies a full event reconstruction. At the hardware trigger stage, events are 
required to have a muon with high px or a hadron, photon or electron with high transverse 
energy in the calorimeters. For hadrons, the transverse energy threshold is 3.5 GeV. The 
software trigger requires a two-, three- or four-track secondary vertex with a signihcant 
displacement from the primary pp interaction vertices (PVs). At least one charged particle 
must have a transverse momentum px > 1.7 GeV/c and be inconsistent with originating 


from a PV. A multivariate algorithm 64 is used for the identihcation of secondary vertices 
consistent with the decay of a 6 hadron. The px of the photon from D*~ decay is too low 
to contribute to the trigger decision. 

Simulated events are used to characterise the detector response to signal and certain 


types of background events. In the simulation, pp collisions are generated using Pythia 65 
with a specihc LHGb conhguration 66 . Decays of hadronic particles are described 
by EvtGen 67 , in which hnal-state radiation is generated using Photos 68 . The 


interaction of the generated particles with the detector, and its response, are implemented 


3 











using the Geant4 toolkit 69 as described in Ref. 70 


3 Event selection 


Signal candidates are formed by combining candidates, reconstructed in the decay 
channel K^tt~ , with two additional pion candidates of opposite charge. Reconstrncted 
tracks are required to be of good quality and to be inconsistent with originating from a PV. 
They are also required to have sufficiently high p and pt and to be within kinematic regions 
where reasonable particle identihcation (PID) performance is achieved, as determined 
by calibration samples of D*~^ —)■ D^(K~n~^)7r~^ decays. The fonr hnal state tracks are 
required to be positively identified by the PID system. The daughters are required to 
form a good quality vertex and to have an invariant mass within 100 MeV/c^ of the known 
mass [^. The candidates and the two charged pion candidates are reqnired to 
form a good vertex. The reconstructed and vertices are required to be significantly 
dis plac ed from the PV. To improve the B^ candidate invariant mass resolution, a kinematic 
is used, constraining the candidate to its known mass 


£t 71 
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By reqniring the reconstrncted vertex to be displaced downstream from the re¬ 
constrncted B^ vertex, backgronnds from both charmless B decays and direct prompt 
charm production coming from the PV are reduced to a negligible level. Background from 
Il*(2010)“ decays is removed by reqniring > 2.1GeV/c^. Backgronnds from 

doubly mis-identified —)■ K^ti~ or doubly Cabibbo-suppressed —)■ K~'k^ decays are 
also removed by this requirement. 

To further distinguish signal from combinatorial background, a multivariate analysis 
based on a Fisher discriminant (7^ is applied. The sPlot techniqne |73 is nsed to 


statistically separate signal and background events with the B^ candidate mass used as the 
discriminating variable. Weights obtained from this procedure are applied to the candidates 
to obtain signal and background distributions that are used to train the discriminant. The 
Fisher discriminant uses information about the event kinematic properties, vertex quality, 
IP and Pt of the tracks and flight distance from the PV. It is optimised by maximising 
the pnrity of the signal events. 

Signal candidates are retained for the Dalitz plot analysis if the invariant mass of the 
B^ meson lies in the range [5250, 5310] MeV/c^ and that of the meson in the range 
[1840, 1890] MeV/c^ (called the signal region). Once all selection requirements are applied, 
less than 1 % of the events contain multiple candidates, and in those cases one candidate 
is chosen randomly. 

Backgronnd contribntions from decays with the same topology, bnt having one or 
two mis-identihed particles, are estimated to be less than 1 % and are not considered in 
the Dalitz analysis. These backgronnd contributions include decays like B^ —, 
B2 D^K- 


K+K 


7T 


74 


—)■ D^pn 


75 


and B^ —)■ 


TT 


with -)■ TT’' 


TT or 




Partially reconstructed decays of the type B^ —)■ X, where one or more particles 

are not reconstructed, have similar efficiencies to the signal channel decays. They are 
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5250 5300 5350 5400 5450 5500 

m(D\+7i“) [MeV/c^] 

Figure 2: Invariant mass distribution of —)• candidates. Data points are shown in 

black. The fit is shown as a solid (red) line with the background component displayed as dashed 
(green) line. 


distributed in the region below the mass. By requiring the invariant mass of 
candidates to be larger than 5250 MeV/c^, these backgrounds are reduced to a negligible 
level, as determined by simulated samples of B^ —)■ and B^ —)■ D*^p{770) with 

D*^ decaying into or D^ti^ under different hypotheses for the D*^ helicity. 

The signal and combinatorial background yields are determined using an unbinned 
extended maximum likelihood ht to the invariant mass distribution of B^ candidates. The 
invariant mass distribution is shown in Fig. with the ht result superimposed. The ht 
uses a Crystal Ball (CB) function 76 convoluted with a Gaussian function for the signal 
distribution and a linear function for the combinatorial background distribution in the 
mass range of [5250, 5500] MeV/c^. Simulated studies validate this choice of signal shape 
and the tail parameters of the CB function are hxed to those determined from simulation. 
Table [T] summarises the ht results on the free parameters, where pso is the mean peak 
position and ctg is the width of the Gaussian function. The parameter ctcb is the width of 
the Gaussian core of the CB function. The parameters /cb and pi give the ht fraction 
of the CB function and the slope of the linear function that describes the background 
distribution. The yields of signal (z/°) and background events given in Table are 
calculated within the signal region. The purity is (97.8 ± 0.2) %. 


4 Dalitz plot analysis formalism 


The analysis of the distribution of decays across the Dalitz plot 15 allows a determination 
of the amplitudes contributing to the three-body B^ —)■ decay. Two of the three 
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Table 1: Results of the fit to the invariant mass distribution of —)• D^'k^-k candidates. 

Uncertainties are statistical only. 


Parameter 


Value 

/i^o 

5282.1 

± 

0.2 MeV/c^ 

0-G 

33.6 

± 

5.4 MeV/c^ 

O'CB 

13.4 

± 

0.3 MeV/c2 

fcB 

0.908 

± 

0.025 

Pi 

-0.152 

± 

0.035 (GeV/c2)-i 


9565 

± 

116 


215 

± 

19 


possible two-body invariant mass-squared combinations, which are connected by 

+ m^(D°7r“) -|- = m\o + m|,o + 2m^, (1) 

are sufficient to describe the kinematics of the system. The two observables 
and where resonances are expected to appear, are chosen in this paper. These 

observables are calculated with the masses of the and mesons constrained to their 
known values The invariant mass resolution has negligible effect and therefore it is 
not modeled in the Dalitz plot analysis. 

The total decay amplitude is described by a coherent sum of amplitudes from resonant 
or nonresonant intermediate processes as 

M{x) = ^CiAi{x). (2) 

i 


The complex coefficient Cj and amplitude Ai{x) describe the relative contribution and 
dynamics of the Tth intermediate state, where x represents the (m^(T)°7r“), m^(7r’''7r“)) 
coordinates in the Dalitz plot. The Dalitz plot analysis determines the coefficients q. In 
addition, ht fractions and interference ht fractions are also calculated to give a convention- 
independent representation of the population of the Dalitz plot. The ht fractions are 
dehned as 


/ \ciAi{x)\‘^ dx 
1 1 J2iCiAi{x)\'^ dx' 


( 3 ) 


and the interference ht fractions between the resonances i and j {i < j) are dehned as 


T - = 

ij 


/2Re[, 


CiC'jAi{x)A* 


[x)]dx 


I I 


( 4 ) 


where the integration is performed over the full Dalitz plot with m(D°7r^) >2.1 GeV/c^. 
Due to these interferences between different contributions, the sum of the ht fractions is 
not necessarily equal to unity. 
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( 5 ) 


The amplitude Ai(x) for a specific resonance r with spin L is written as 
Ai{x) = Fi^^\q,qo) x X Tl{x) x R{x). 


The functions F^\q,qQ) and Fr^\p,po) are the Blatt-Weisskopf barrier factors 


AL), 


77 


for 


the production, —)■ rh^, and the decay, r —)■ hih 2 , of the resonance, respectively. The 
parameters p and q are the momenta of one of the resonance daughters (hi or ^ 2 ) and of 
the bachelor particle (h^), respectively, both evaluated in the rest frame of the resonance. 
The value po {qo) represents the value of p {q) when the invariant mass of the resonance is 
equal to its pole mass. The spin-dependent Fb and Fr functions are defined as 


L = 0 
L = 1 

L = 2 

L = 3 

L = A 


: F(o)(z,Zo) = l, 


1 + ^ ’ 



{^0 ~ 3 )^ + 9^0 
{z-3y + 9z ’ 


( 6 ) 


where Z{q) is equal to (rew x q'(o))^ or (rsw x P{p)Y 
resonance, rewi is taken to be 1.6 GeV“^ x he (= 0.3 fm) [78 


;^o(;^o-15)2 + 9(2^0-5) 
z{z - 15y + 9{2z - 5) ’ 

{z^ - 45zo + 105)2 + 25zo{2zo - 21)2 
(^2 - 45z + 105)2 + 25z(2z - 21)2 ’ 

The value for the radius of the 


The function Tb{x) represents the angular distribution for the decay of a spin L 
resonance. It is dehned as 


L = 0 
L = 1 
L = 2 
L = 3 
L = 4 


^0 — 1, 

Ti = \/l + cos 9 X qp, 

T2 = (|/2 + 3/2)(cos2 0-l/3)xgV, 

T 3 = a/ 1 + ?/2(l + 2|/^/5)(cos^ 9 — 3 cos(9)/5) x q^p^, 

T 4 = ( 81 /V 35 + AOy^/35 + l)(cos^0 - 3Ocos2(0)/35 + 3/35) x gV 


(7) 


4„4 


The helicity angle, 6, of the resonance is defined as the angle between the direction of the 
momenta p and q. The y dependence accounts for relativistic transformations between the 
and the resonance rest frames 79,80 , where 


^ 2 mAa + rn?{hih2) 


m 


hs 


(g) 

2m{hih2)mB0 

Finally, R{x) is the resonant lineshape and is described by the relativistic Breit-Wigner 
(RBW) function unless specihed otherwise, 

1 


RBW(s) = 


ml — s — imj.rR)(s) ’ 


(9) 
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where s = m?‘{hih 2 ) and is the pole mass of the resonance; the mass-dependent 

width, is dehned as 




( 10 ) 


where Fq is the partial width of the resonance, i.e., the width at the peak mass s = rrir. 
The lineshapes of p(770), p(1450) and p(1700) are described by the Gounaris-Sakurai 


(GS) function 81 


GS(s) = 


ml{l + T^g/mr) 
ml — s + f{s) — imrTp{s) ’ 


( 11 ) 


where 


fis) 

h{s) 


ml 
J- 0^ 
Po 


{h{s) - h{ml))p‘^ + {ml - 


9 = 


2 p 

-^log 

vr a/s 

3 ml 

-^log 

TT Po 



3 

r 9 “1 

[po\ 


s 


\/s + 2p 
2m,r 

mr + 2po 

2m^ 

1/2 


+ 


mr 


2'kpo 


mlmr 

m 


and rp(s) = Fc 

The p — uj interference is taken into account by 

Rp-uj{s) = GSp(77o)(s) X (1 -f- ae®®RBW^(782)(s)), 


( 12 ) 


(13) 


where Fq is used, instead of the mass-dependent width F^^/s), for ci;(782) 82 


The 77* (2010) contribution is vetoed as described in Sec. Possible remaining 
contributions from the Zl*(2010)“ RBW tail or general P-waves are modelled as 




Rii*(20i0)("i 7r+)) = e 


(14) 


where /3i and (32 are free parameters. 

The 7r'''7r“ S-wave contribution is modelled using two alternative approaches, the 
Isobar model 16-18 or the K-matrix model 1^. Gontributions from the /o(500), /o(980). 


/o(2020) resonances and a nonresonant component are parametrised separately in the 
Isobar model and globally by one amplitude in the K-matrix model. 

In the Isobar model, the /o(2020) resonance is modelled by a RBW function and the 
modelling of the /o(500), /o(980) resonances and the nonresonant contribution are described 
as follows. The Bugg resonant lineshape 83 is employed for the /o(500) contribution, 

(15) 


^/o(50l))(^) — 1/ 


9l 2 
mi 


sa 


sa 


z{s) - imrTtot{s) 
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where 


mrTi{s) 

9 iis) 


= 9i- 


sa 


Pi{s), 


Z(S = 


mr{bi + b 2 s) exp(-(s - ml)/A), 

3 i{s) 


jii.'ml), 


Ji(s) = - 


2 + Pi log 


1 -pi 


1 + Pi 

mrT 2 {s) = 0.6pi(s)(s/m^) exp(-a|s - 4 m|-|)p 2 (s), 
mrTsls) = 0.2gl{s){s/ml)exp{-a\s-AmlDpsis), 

mrT^^s) = mrg4nP4n{s)/p4n{ml), 

P 4 ^(s) = 1/[1 + exp(7.082 - 2.845s)], 


and rtot(s) = y^r^(s)- 


2 = 1 


(16) 


The parameters are hxed to = 0.953 GeV/c^, s^ = 0.41 bi = 1.302 GeVyc"^, 
62 = 0.340, A = 2.426 GeVyc"^ and = 0.011 GeV/c^ j^. The phase-space factors of the 
decay channels mi, KK and pp correspond to pi, 2 , 3 (s), respectively and are dehned as 


y / 

Pl,2,3('S) = 1/1—4- 


m 


1,2,3 


1, 2, and 3 = 7r,K and p. 


(17) 


The Flatte formnla 84 is nsed to describe the /o(980) lineshape, 

1 


where 


Rfo{980){s) — 


ml - s- imr(p 7 r 7 r(s)pi + PKK{s)g 2 ) ’ 


p7T7v{^) 3 ^ ^ 




4m£+ 1 L 4m^o 




(18) 


(19) 


The parameters pi( 2 ) and 46 are = 939.9 ± 6.3MeV/c^, pi = 199 ± 30MeV and 
5'2/fi'i = 3.0 ± 0.3. 

The nonresonant contribntion is described by 


RNR{m\7r+7r-),m\D%+)) = ). 


( 20 ) 


Its modnlns eqnals nnity, and a slowly varying phase over m?{Ti^Ti ) acconnts for rescat¬ 
tering effects of the tt+tt" hnal state and a is a free parameter of the model. 
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The K-matrix formalism describes the production, rescattering and decay of the 
7 r''' 7 r“ S-wave in a coherent way. The scattering matrix S, from an initial state to a hnal 
state, is 

S = I + 2z{p^Y^^Tp^/^, ( 21 ) 

where I is the identity matrix, p is a diagonal phase-space matrix and T is the transition 
matrix. The unitarity requirement S'S'l = I gives 

(T-^ + tpY = T-^ + ip. (22) 

The K-matrix is a Lorentz-invariant Hermitian matrix, dehned as = T~^ -|- ip. The 
amplitude for a decay process, 

A. = (I- iKp)pPi, (23) 

is computed by combining the K-matrix obtained from scattering experiment with a 
production vector to describe process-dependent contributions. The K-matrix is modelled 
as a hve-pole structure. 


K«(s) 



g.°9° 

— s 


1 ^scatt 

I rscatt ^ ^0 

“r J ij 


S — S 


scatt 
0 


1 - 5A0 / _ Sa^'\ 

s - sao V 2 / ’ 


(24) 


where the indexes i,j = 1,2, 3,4,5 correspond to hve decay channels: vrTr, KK, pp, pp' 
and multi-meson (mainly dvr states) respectively. The coupling constant of the bare state 
a to the decay channel i, gf, is obtained from a global fit of scattering data and is listed in 
Table The mass rria is the bare pole mass and is in general different from the resonant 
mass of the RBW function. The parameters /fj^*^** and are used to describe smooth 


s — 


85,86 


scattering processes. The last factor of the K-matrix, 

singularities near the tt+tt" threshold, the so-called 
property of the K-matrix imposes the relation j 

are considered, \i i ^ 1 and j 7 ^ 1, is set to 0. The production vector is modelled 
with 


1 — 

S—.SAO 

“Adler zero” 

pscatt _ Mscatt 

J i 


SAm^ 


, regulates the 

The Hermitian 
and since only tt+tt” decays 




1 _ „prod p a 

^prod *0 \ ^ ycuUj 

Pi „ 2-^ yWl _ g 

S - Sq ^ lll'a * 


(25) 


where and 13a are free parameters. The singularities in the K-matrix and the 

production vector cancel when calculating the amplitude matrix element. 


5 Dalitz plot fit 


An unbinned extended maximum likelihood fit is performed to the Dalitz plot distribution. 
The likelihood function is dehned by 


C = C' X 


A/2TTas 


exp 


- zz, 


0'i2 


2 a2 


X 


V^(Jb 


exp 


2al 


(26) 
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Table 2: The K-matrix parameters used in this paper are taken from a global analysis of 7r'*“7r 
scattering data 22 . Masses and coupling constants are in units of GeV/c^. 


rUa 

9"+^- 

„a _ 
^KK 

^47r 

JTjr) 

^7)7]' 

0.65100 

0.22889 

-0.55377 

0.00000 

-0.39899 

-0.34639 

1.20360 

0.94128 

0.55095 

0.00000 

0.39065 

0.31503 

1.55817 

0.36856 

0.23888 

0.55639 

0.18340 

0.18681 

1.21000 

0.33650 

0.40907 

0.85679 

0.19906 

-0.00984 

1.82206 

0.18171 

-0.17558 

-0.79658 

-0.00355 

0.22358 


fscatt 

Jll 

/scatt 

Jl2 

/scatt 

il3 

.f scatt 
Jl4 

/scatt 

Jib 


0.23399 

0.15044 

-0.20545 

0.32825 

0.35412 


c>scatt 

-^0 

prod 

-^0 


sa 



-3.92637 

-3.0 

-0.15 

1 



where 


a = 


iv! 


N r 


n 

2 = 1 


Vs + 


fs T 


Vb 


Vs + Vb 


9b) 


(27) 


The background probability density function (PDF) is given by /^(x; 6b) and is described 


in Sec. 5.1 The signal PDF, fs{xi',6s), is described by 

M{xi;9s)e{xi) 
f M(x; 6s)e{x)dx' 


(28) 


where the decay amplitude, M(x;9s), is described in Sec. |^and the efficiency variation 
over the Dalitz plot, e{x), is described in Sec. 5.2 The £t parameters, 9s and 9b, include 
complex coefficients and resonant parameters like masses and widths. The valne N is 
the total nnmber of reconstructed candidates in the signal region. The number of signal 
and backgronnd events, Ug and Ub, are floated and constrained by the yields, and 
determined by the D^7r~^7i~ mass ht and shown with their statistical nncertainties in 
Tabled 


5.1 Background modelling 

The only signihcant source of candidates in the signal region, other than 
decays, is from combinatorial backgronnd. It is modelled nsing candidates in the npper 
m{D^7T~^7r~) sideband ([5350, 5450] MeV/c^) with a looser reqnirement on the Fisher 
discriminant, and is shown in Fig. The looser requirement gives a similar distribution in 
the Dalitz plane but with lower statistical fluctuations. The Dalitz plot distribution of the 
combinatorial background events lying in the upper-mass sideband is considered to provide 
a reliable description of that in the signal region, as no dependence on m{D^7r^7T~) is 
found by studying the Dalitz distribution in a different upper-mass sideband region. The 


combinatorial background is modelled with an interpolated non-parametric PDF 87 
using an adaptive kernel-estimation algorithm (Ml. 
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Figure 3: Density profile of the combinatorial background events in the Dalitz plane obtained 
from the upper m{D^7r~^7r~) sideband with a looser selection applied on the Fisher discriminant. 


5.2 Efficiency modelling 


The efficiency function e{x) accounts for effects of reconstruction, triggering and selection 
of the —)■ D^7r~^7T~ signal events, and varies across the Dalitz plane. Two simulated 

samples are generated to describe its variation with several data-driven corrections. One 
is uniformly distributed over the phase space of the Dalitz plot and the other is uniformly 
distributed over the square Dalitz plot, which models efficiencies more precisely at the 
kinematic boundaries. The square Dalitz plot is parametrised by two variables m' and 6' 
that each varies between 0 and 1 and are defined as 


, i m[7l '71 - m{7l ' 71 Jmin 

m = — arccos /——r -—-—r— 

TT V "ifTT+TT )max “ m['K+'K jmin 


where m(7r’''7r“)max = ~ and 6{7i~^7r~) is the helicity angle 

of the TT+TT" system. 

The two samples are fitted simultaneously with common fit parameters. A 4th-order 
polynomial function is used to describe the efficiency variation over the Dalitz plot. As the 
efficiency in the simulation is approximately symmetric over m^(D° 7 r+) and m^(D° 7 r“), 
the polynomial function is defined as 


e{x,y) oc + aQ{x + y) + ai{x + yY + a 2 {xy) + a'i{x + yf 

- 1 - 04 ( 2 ; -|- y)xy + a^{x -f yY + + yYxy + ajx'^y'^, (30) 


where 


m^(D°7r+) — ml 
{mso — — ml 


and 


y = 


m?{D^77 ) — ml 
{mso — m^^y — ml ’ 


(31) 


with ml defined as [{m£,o -fmjr)^ -t- (m^o — mT^y]/2. The fitted efficiency distribution over 
the Dalitz plane is shown in Fig. 
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Figure 4: Efficiency function for the Dalitz variables obtained in a fit to the LHCb simulated 
samples. 


The efficiency is corrected using dedicated control samples with data-driven methods. 
The corrections applied to the simulated samples include known differences between 
simulation and data that originate from the trigger, PID and tracking. 

5.3 Results of the Dalitz plot fit 

The Dalitz plot distribution from data in the signal region is shown in Fig. The analysis 
is performed using the Isobar model and the K-matrix model. The nominal ht model in 
each case is defined by considering many possible resonances and removing those that do 
not significantly contribute to the Dalitz plot analysis. The resulting resonant contributions 
are given in Table while the projections of the fit results are shown in Fig. (Fig. for 
the Isobar (K-matrix) model. 

The comparisons of the S-wave results for the Isobar model and the K-matrix model are 
shown in Fig. The results from the two models agree reasonably well for the amplitudes 
and phases. In the tt+tt” mass-squared region of [1.5,4.0] GeV^/c"^, small structures are 
seen in the K-matrix model, indicating possible contributions from /o(1370) and /o(1500) 
states. These contributions are not signihcant in the Isobar model and are thus not 
included in the nominal ht: adding them results in marginal changes and shows similar 
qualitative behaviour to the K-matrix model as displayed on Fig. The measured S-waves 
from both models qualitatively agree with predictions given in Ref. [Mj . 

To see more clearly the resonant contributions in the region of the p(770) resonance, 
the data are plotted in the tt+tt” invariant mass-squared region [0.0, 2.1] GeV^/c^ in Fig. 

In the region around 0.6 GeV^/c^, interference between the p(770) and Ci;(782) resonances 
is evident. In the tt+tt” S-wave distributions of both the Isobar model and the K-matrix 
model, a peaking structure is seen in the region [0.9,1.0] GeV^/c"^, which corresponds to 
the /o(980) resonance. The structure in the region [1.3,1.8] GeV^/c^ corresponds to the 
spin-2 /2(1270) resonance. 
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Figure 5: Dalitz plot distribution of candidates in the signal region, including background 
contributions. The red line shows the Dalitz plot kinematic boundary. 


Table 3: Resonant contributions to the nominal fit models and their properties. Parameters and 


uncertainties of p{770), w(782), /9(1450) and /9(1700) come from Ref. 90 , and those of /2(1270) 
and /o(2020) come from Ref. [32] . Parameters of /o(500), /o(980) and K-matrix formalism are 
described in Sec. 


Resonance 


Spin Model (MeV/c^) Fq (MeV) 


p(770) 

u{782) 

p(1450) 

p(1700) 

/2(1270) 


TTTT 


-wave 


/o(500) 

/o(980) 

/o(2020) 

Nonresonant 


D^ti P-wave 

1 

Eq. 14 

Floated 

±>*(2400)- 

0 

RBW 

Floated 

D*(2460)- 

2 

RBW 

Floated 

D}(2760)- 

3 

RBW 

Floated 


1 

1 

1 

1 

2 


GS 
Eq, 
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GS 

GS 

RBW 


775.02 ± 0.35 
781.91 ±0.24 
1493 ± 15 
1861 ± 17 
1275.1 ± 1.2 


149.59 ±0.67 
8.13 ±0.45 
427 ± 31 
316 ± 26 
185.1 


0 K-matrix 


See Sec. 


0 

0 

0 

0 


Eq. ^ 
Eq. 18 
RBW 


Eq. 20 


See Sec. 4 
See Sec. 4 
1992 ± 16 442 ± 60 

See Sec. El 


2.9 

2.4 


Distributions in the invariant mass-squared region [6.4,10.4] GeV^/c^ of 
are shown in Eig. 10 There is a signihcant contribution from the D}(2760)“ resonance 
observed in Ref. and a spin-3 assignment gives the best description. A detailed 


discussion on the determination of the spin of Zl}(2760) is provided in Sec. 8.2 

The ht quality is evaluated by determining a value by comparing the data and the 
ht model in Wins = 256 bins that are dehned adaptively to ensure approximately equal 
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Figure 6: Projections of the data and Isobar fit onto (a) m^(7r'*“7r“) and (c) m?‘{D^'K~) with a 
linear scale. Same projections shown in (b) and (d) with a logarithmic scale. Components are 
described in the legend. The lines denoted D^tt~ and 7r'’'7r“ include the coherent sums of all 
D^n- resonances, 7r'*“7r resonances, and tt+tt S-wave resonances. The various contributions do 
not add linearly due to interference effects. 


population with a minimum bin content of 37 entries. A value of 287 (296) is found for the 
Isobar (K-Matrix) model based on statistical uncertainties only. The effective number of 
degrees of freedom (nDoF) of the is bounded by A^bins — 1 and A^bins — Alpars — 1, where 
Apars is the number of parameters determined by the data. Pseudo experiments give an 
effective number of 234 (235) nDoF. 

Further checks of the consistency between the fitted models and the data are performed 
with the unnormalised Legendre polynomial weighted moments as a function of m?{D^Ti~) 
and m?{TT~^TT~). The corresponding distributions are shown in Appendix [A} 
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Figure 7: Projections of the data and K-matrix fit onto (a) m^(7r+7r“) and (c) m‘^(D^7r~) with a 
linear scale. Same projections shown in (b) and (d) with a logarithmic scale. Components are 
described in the legend. The lines denoted D^tt~ and 7r'’'7r“ include the coherent sums of all 
D^n- resonances, 7r'*“7r resonances, and tt+tt S-wave resonances. The various contributions do 
not add linearly due to interference effects. 


6 


Measurement of the 
tion 


branching frac- 


Measuring the branching fractions of the different resonant contributions requires knowledge 
of the branching fraction. This branching fraction is normalised relative to 

the —)■ Zl)*(2010)“7r’^ decay that has the same final state, so systematic uncertainties are 
reduced. Identical selections are applied to select B^ —)■ Zl)*(2010)“7r’'' and B^ —)■ D^tt^ 


IT 


candidates, the only difference being that m{D^'K~) < 2.1 GeV/c^ is used to select 
D*{2010)~ candidates. The kinematic constraints remove backgrounds from doubly 
mis-identified —)■ K^'k~ or doubly Cabibbo-suppressed —)■ K~7i^ decays and no 
requirement is applied on 

The invariant mass distributions of m(D^7r~) and m(D^n~^7r~) for the B^ —)■ 
T)*(2010)“7r+ candidates are shown in Fig. 11 and are fitted simultaneously to deter¬ 
mine the signal and background contributions. The T)*(2010)“ signal distribution is 
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Figure 8: Comparison of the 7r'’'7r“ S-wave obtained from the Isobar model and the K-matrix 
model, for (a) amplitudes and (b) phases. The K-matrix model is shown by the red solid line, 
two scenarios for the Isobar model with (black long dashed line) and without (blue dashed line) 
/o(1370) and /o(1500) are shown. 




Figure 9: Distributions of m^(7r'*“7r“) in the /9(770) mass region. The different fit components 
are described in the legend. Results from (a) the Isobar model and (b) the K-matrix model are 
shown. 


modelled by three Gaussian functions to account for resolution effects while its background 
is modelled by a phase-space factor. The modelling of the signal and background shapes 
in the m{D^TT^'K~) distribution are described in Sec. The —)■ T)*(2010) 7r+ yield in 
the signal region is 7327 ± 85. 

The efficiencies for selecting T)*(20I0)“7r’'' and B^ —)■ decays are 

obtained from simulated samples. To take into account the resonant distributions in the 
Dalitz plot, the B^ —)■ simulated sample is weighted using the model described in 

the previous sections. The average efficiencies are (1.72±0.05) x 10“^ and (4.96±0.05) x 10“"^ 
for the —)■ T)*(2010)“7r+ and B^ —)■ decays. 

Using the branching fractions of i3(i?° —)• T)*(20I0)“7r+) = (2.76 ± 0.13) x 10“^ 
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Figure 10: Distributions of rin?{D^'K~) in the Ilj(2760)“ mass region. The different fit components 
are described in the legend. Both results from (a) the Isobar model and (b) the K-matrix model 
are shown. 



Figure 11: Invariant mass distributions of (a) and (b) m{D^TT~) for —>• 

candidates. The data is shown as black points with the fit superimposed as 

red solid lines. 


and i3(L)*(2010) 




= (67.7 ± 0.5)% 


32 


the derived branching fraction of 
TT in the kinematic region m{D^TT^) >2.1 GeV/c^ is (8.46±0.14±0.40) x 10“^, 


where the hrst nncertainty is statistical and the second uncertainty comes from the 
branching fraction of the normalisation channel. 


7 Systematic uncertainties 

7.1 Common systematic uncertainties and checks 

Two categories of systematic uncertainties are considered, each of which is quoted separately. 
They originate from the imperfect knowledge of the experimental conditions and from the 
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assumptions made in the Dalitz plot fit model. The Dalitz model-dependent uncertainties 
also account for the precision on the external parameters. The various sources are assumed 
to be independent and summed in quadrature to give the total. 

Experimental systematic uncertainties arise from the efficiency and background mod¬ 
elling and from the veto on the T)*(2010)“ resonance. Those corresponding to the signal 
efficiency are due to imperfect estimations of PID, trigger, tracking reconstruction effects, 
and to the hnite size of the simulated samples. Each of these effects is evaluated by the 
differences between the results using efficiencies computed from the simulation and from 
the data-driven methods. The systematic uncertainties corresponding to the modelling of 
the small residual background are estimated by using different sub-samples of backgrounds. 
The systematic uncertainty dne to the veto on the T)*(2010)“ resonance is assigned by 
changing the selection reqnirement from m‘^{D^7T^) > 2.10GeV/c^ to 2.05GeV/c^. 

The systematic uncertainties related to the Dalitz models considered (see Sec.|^ include 
effects from other possible resonant contributions that are not included in the nominal ht, 
from the modelling of resonant lineshapes and from imperfect knowledge of the parameters 
of the modelling, i.e., the masses and widths of the 7r’''7r“ resonances considered, and the 
resonant radius. 

The non-signfficant resonances added to the model for systematic stndies are the 
/o(1300), /o(1500), /2(1525), and D*(2650)“ (/2(1525) and D*(2650)“) mesons for the 


Isobar (K-matrix) model (^|^|^|^. The spin of the D*(2650)“ resonance is set to 1. 
The differences between each alternative model and the nominal model are conservatively 
assigned as systematic nncertainties. 

The radins of the resonances (rew) is set to a nniqne valne of 1.6GeV“^ x he in the 
nominal ht. In the systematic stndies, it is hoated as a free parameter and its best ht 
value is 1.84 ± 0.05 GeV“^ x he (1.92 ± 0.31 GeV“^ x he) for the Isobar (K-matrix) model. 
The value 1.85 GeV“^ x he is chosen to estimate the systematic nncertainties dne to the 
imperfect knowledge of this parameter. 

The masses and widths of the tt+tt” resonances considered are treated as free parameters 
with Ganssian constraints according to the inpnts listed in Table The diherences 
between the results from those hts and those of the nominal hts are assigned as systematic 
uncertainties. 

Eor the Isobar model, additional systematic nncertainties dne to the modelling of 
the /o(500) and /o(980) resonances are considered. The Bngg model 83 for the /o(500) 
resonance and the Elatte model 


84 


for the /o(980) resonance, nsed in the nominal ht, 
are replaced by more conventional RBW fnnetions. The masses and widths, left as 
free parameters, give 553 ± 15MeV/c^ and 562 ± 39MeV, for the /o(500) meson and 
981 ± 13MeV/c^ and 191 ± 39MeV, for the /o(980) meson. The resnlting diherences to 
the nominal ht are assigned as systematic nncertainties. 

The kinematic variables are calcnlated with the masses of the and mesons 


constrained to their known values 32 . These kinematic constraints ahect the extraction of 


the masses and widths of the resonances. The current world average value for the 

meson mass is 5279.58 ± 0.17MeV/c^ and for the meson is 1864.84 ± 0.07MeV/c^ 32 


A conservative and direct estimation of the systematic nncertainties on the masses and 
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Table 4: Systematic uncertainties on B{B^ —)• D^-k^'k ). 


Source 

Uncertainty (xlO 

PID 

0.02 

Trigger 

0.13 

Reconstruction 

< 0.01 

Size of simulated sample 

0.26 

D*(2010)“ mass model 

< 0.01 

Dalitz structure 

0.04 

Total 

0.29 


widths of the D^tt~ resonances is provided by the snm in qnadratnre of the and 
mass nncertainties. The effects of mass constraints are fonnd to be negligible for the ht 
fractions, modnli and phases of the complex coefficients. 

The systematic uncertainties are summarised for the Isobar (K-matrix) model Dalitz 
analysis in Appendix Systematic uncertainties related to the measurements performed 
with the Isobar formalism are listed in Tables [T^ to [TTI while those for the K-matrix 


formalism are given in Tables to [21 
uncertainties are due to the T)*(2010) 


In most of cases, the dominant systematic 


veto and the model uncertainties related to other 
resonances not considered in the nominal ht. In the Isobar model, the modelling of the 
/o(500) and /o(980) resonances also have non-negligible systematic effects. 

Several cross-checks have been performed to study the stability of the results. The 
analysis was repeated for different Fisher discriminant selection criteria, different trigger 
requirements and different sub-samples, corresponding to the two data-taking periods and 
to the two half-parts of the invariant mass signal region, above and below the 

Results from those checks demonstrate good consistency with respect to the 


mass 


32 


nominal ht results. No bias is seen, therefore no correction is applied, nor is any related 
uncertainty assigned. 


7.2 Systematic uncertainties on the D^7r‘'‘7r branching 

fraction 

The systematic uncertainties related to the measurement of the —)■ branching 

fraction are listed in Table [^ The systematic uncertainties on the PID, trigger, reconstruc¬ 
tion and statistics of the simulated samples are calculated in a similar way to those of the 
Dalitz plot analysis. Other systematic uncertainties are discussed below. 

The systematic uncertainty on the modelling of the D^tt~ and D*(2010)“7r’'' invariant 
mass distributions is estimated by counting the number of signal events in the B^ signal 
region assuming a hat background contribution. The D*(2010)“ mass region is restricted 
to the range [2007, 2013] MeV/c^ for this estimate. The calculated branching fraction is 
nearly identical to that from the mass ht and thus has a negligible contribution to the 
systematic uncertainty. The signal purity of B^ —)■ D*(2010)“7r’'' is more than 99%. 
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Table 5: Statistical significance (a) of 7r'*“7r“ resonances in the Dalitz plot analysis. For the 
statistically significant resonances, the effect of adding dominant systematic uncertainties is 
shown (see text). 


Resonances 

w(782) 

/o(980) 

/o(1370) 

p(1450) 

/o(1500) 

/'(1525) 

p(1700) 

/o(2020) 

Isobar 

8.0 

10.7 

1.1 

8.7 

1.1 

3.6 

4.5 

10.2 

K-matrix 

8.1 

n/a 

n/a 

8.6 

n/a 

2.6 

2.2 

n/a 

With syst. 

7.7 

7.0 

n/a 

8.7 

n/a 

n/a 

n/a 

4.3 


To account for the effect of resonant structures on the signal efficiency, the data sample 
is divided using an adaptive binning scheme. The average efficiency is calculated in a 
model independent way as 

V N- 




where Ni is the number of events in bin i and e* is the average efficiency in bin i calcnlated 
from the efficiency model. The difference between this model-independent method and 
the nominal is assigned as a systematic uncertainty. 


8 Results 


8.1 Significance of resonances 

The Isobar and K-matrix models employed to describe the Dalitz plot of the D^7T~^'n'~ 

decay include all of the resonances listed in Table The statistical significances of well- 
established 7r’''7r“ resonances are calculated directly with their masses and widths hxed 
to the world averages. They are computed as the relative change of the minimum of the 
negative logarithm of the likelihood (NLL) function with and without a given resonance. 
Besides the tt+tt" resonances listed in Table the significances of the /o(1370), /o(1500) 
and /2(1525) are also given. The results, expressed as multiples of Gaussian standard 
deviations (cr), are summarised in Table All of the other tt+tt” resonances not listed in 
this Table have large statistical significances, well above fivestandard deviations. 

To test the signihcance of the D}(2760)“ state, where J = 3 (see Sec. 8.2), an ensemble 


reported in Ref. 29 . The two other orbitally 


of pseudo experiments is generated with the same number of events as in the data sample, 
using parameters obtained from the ht with the D}(2760)“ resonance excluded. The 
difference of the minima of the NLL when htting with and without D}(2760)“ is used as 
a test statistic. It corresponds to 11.4 a (11.5 a) for the for the Isobar (K-matrix) model 
and confirms the observation of D 
excited D resonances, D}(2650) 
the same paper, are added into the nominal fit model with different spin hypotheses and 
tiny improvements are found. They also do not describe the data in the absence of the 
Dj{2760)~. Those resonances are thus not confirmed by this analysis. Finally, an extra 
D^'k~ resonance, with different spin hypotheses (J = 0,1, 2, 3,4) and with its mass and 


,j(2760) 
and D 


j(3000) , whose observations are presented in 
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width allowed to vary, is added to the nominal ht model and no signihcant contribution is 
found. 

The signihcance of each of the signihcant a;(782), /o(980), p(1450), /o(2020) and 
Dj{2760)~ states is checked while including the dominant systematic uncertainties (see 


Sec. 7.1), namely, the modelling of the /o(500) and /o(980) resonances, the addition of 
other resonant contributions and the modihcation of the 77* (2010)“ veto criteria. In all 
conhgurations, the signihcances of the Ci;(782), /o(980), p(1450) and Z7}(2760)“ resonances 
are greater than 7.7 a, 7.0 a, 8.7 a, and 10.8 a, respectively. The signihcance of the /o(2020) 
drops to 4.3 cr when using a RBW lineshape for the /o(500) resonance. The abundant 
/o(500) contribution is highly signihcant under all of the applied changes. 


8.2 Spin of resonances 


As described in Sec. 5.3, a spin-3 i7}(2760)“ contribution gives the best description of 
the data. To obtain the signihcance of the spin-3 hypothesis with respect to other spin 
hypotheses (J = 0, 1, 2, 4), test statistics are built. Their computations are based on 
the shift of the minimum of the NLL with respect to the nominal ht model when using 
diherent spin hypotheses. The mass and width of the 77}(2760)“ resonance are hoated 
in all the cases. Pseudo experiments are generated using the ht parameters obtained 
using the other spin hypotheses. Signihcances are calculated according to the distributions 
obtained from the pseudo experiments of the test statistic and its values from data. These 
studies indicate that data are inconsistent with other spin hypotheses by more than 10 a. 
Following the discovery of the 77*j (2860)“ meson, which is interpreted as the superposition 


of two particles with spin 1 and spin 3 j^,26 , a similar conhguration for the 77}(2760)“ 
has been tested and is found to give no signihcant improvement in the description of the 
data. To illustrate the preference of the spin-3 hypothesis, the cosine of the helicity angle 
distributions in the mass-squared region of [7.4, 8.2] GeV^/c^ for m^(77°7r“) are shown in 
Fig. 12 under the various scenarios. Based on our result, 77}(2760)“ is interpreted as the 
77}(2760)“ meson. Recently, LHCb observed a neutral spin-1 77*(2760)° state 92 . The 
current analysis does not preclude a charged spin-1 77* state at around the same mass, 
but it is not sensitive to it with the current data sample size. 

Studies have also been performed to validate the spin-0 hypothesis of the 77} (2400)“ 
resonance, as the spin of this state has never previously been conhrmed in experiment |^ . 
When moving to other spin hypotheses, the minimum of the NLL increases by more than 
250 units in all cases, which confirms the expectation of spin 0 unambiguously. 


8.3 Results of the Dalitz plot analysis 

The shape parameters of the 7r’''7r“ resonances are hxed from previous measurements except 
for the nonresonant contribution in the Isobar model. The htted value of the parameter a 


defined in Eq. (20) is —0.363 ± 0.027, which corresponds to a 10 a statistical significance 
compared to the case where there is no varying phase. An expansion of the model by 
including a varying phase in the 77°7r“ axis is also investigated but no signihcantly varying 
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Figure 12: Cosine of the helicity angle distributions in the m?{D^7r~) range [7.4, 8.2] GeV^/c^ 
for (a) the Isobar model and (b) the K-matrix model. The data are shown as black points. The 
helicity angle distributions of the Dalitz plot ht results, without the Ilj(2760)“ and with the 
different spin hypotheses of D}(2760)“, are superimposed. 


phase in that system is seen. The results indicate a weak, but non-negligible, rescattering 
effect in the tt+tt" states, while the rescattering in the D^tt~ states is not signihcant. The 
masses, widths and other shape parameters of the D^7i~ contributions are allowed to 
vary in the analysis. The values of the shape parameters of the D^7t~ P-wave component, 
defined in Eq. Q, are (3i = 0.95 ± 0.05 (0.90 ± 0.04) and /?2 = 0.51 ± 0.06 (0.43 ± 0.05) 
for the Isobar (K-matrix) model. 

The measurements of the masses and widths of the three resonances Z1 q( 2400)“, 
112(2460)“ and 113(2760)“ are listed in Table The present precision on the mass and 
width of the Zlg(2400)“ resonance is improved with respect to Refs. 29,32 
the width of the Zl2(2460)' 


The result for 


meson is consistent with previous measurements, whereas the 
result for the mass is above the world average which is dominated by the measurement using 
inclusive production by LHCb (^. In the previous LHCb inclusive analysis, the broad 
Ilo(2400)“ component was excluded from the ht model due to a high correlation with the 
background lineshape parameters, while here it is included. The present result supersedes 
the former measurement. The Dalitz plot analysis used in this paper ensures that the 
background under the 112(2460)“ peak and the effect on the efficiency are under control, 
resulting in much lower systematic uncertainties compared to the inclusive approach. 

The moduli and the phases of the complex coefficients of the resonant contributions, 
dehned in Eq. ([^, are displayed in Tables and Compatible results are obtained using 
both the Isobar and K-matrix models. The results for the ht fractions are given in Table 
while results for the interference ht fractions are given in Appendix Pseudo experiments 
are used to validate the htting procedure and no biases are found in the determination of 
parameter values. 
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Table 6: Measured masses (m in MeV/c^) and widths (T in MeV) of the -Dq( 2400)“, 1^2(2460)“ 
and 243 ( 2760 )“ resonances, where the hrst uncertainty is statistical, the second and the third are 
experimental and model-dependent systematic uncertainties, respectively. 




Isobar 

K-matrix 


24^(2400) 

m 

2349 ± 6 ± 1 ± 4 

2354 ± 7± 11± 

2 


r 

217± 13± 5± 12 

230 ± 15 ± 18 ± 

11 

24*(2460) 

m 

2468.6 ±0.6 ±0.0 ±0.3 

2468.1 ±0.6 ±0.4 ±0.3 


r 

47.3 ±1.5 ±0.3 ±0.6 

46.0 ±1.4 ±1.7 ±0.4 

24*(2760) 

m 

2798 ± 7± 1± 7 

2802± 11 ± 10 ± 

3 


r 

105± 18± 6± 23 

154 ± 27 ± 13 ± 

9 


Table 7: The moduli of the complex coefficients of the resonant contributions for the Isobar 
model and the K-matrix model. The hrst uncertainty is statistical, the second and the third are 
experimental and model-dependent systematic uncertainties, respectively. 


Resonance 


Isobar ( c* ) 


K-matrix ( Q ) 


Nonresonance 

3.43 ± 

0.22 ± 0.04 ± 

0.51 


n/a 


/ o ( 500 ) 

18.7 ± 

0.70 ± 0.29 ± 

0.80 


n/a 


/ o ( 980 ) 

2.62 ± 

0.25 ± 0.09 ± 

0.46 


n/a 


/ o ( 2020 ) 

4.41 ± 

0.51 ± 0.21 ± 

1.78 


n/a 


p ( 770 ) 


1.0 (fixed) 



1.0 (fixed) 


a ;( 782 ) 

0.30 ± 

0.04 ± 0.00 ± 

0.01 

0.31 ± 

0.04 ± 0.01 ± 

0.01 

/ 9 ( 1450 ) 

0.23 ± 

0.03 ± 0.01 ± 

0.02 

0.28 ± 

0.03 ± 0.08 ± 

0.01 

p ( 1700 ) 

0.078 ± 

0.016 ± 0.006 ± 

0.008 

0.136 ± 0.020 ± 0.077 ± 

0.011 

/ 2 ( 1270 ) 

0.072 ± 

0.002 ± 0.000 ± 

0.005 

0.073 ± 0.002 ± 0.006 ± 

0.003 

24^7r“ P-wave 

18.8 ± 

0.7 ± 0.3 ± 

1.9 

19.6 ± 

0.7 ± 0.7 ± 

0.6 

24 ( 5 ( 2400 )“ 

12.1 ± 

0.8 ± 0.3 ± 

0.6 

13.1 ± 

1.0 ± 0.8 ± 

0.5 

24 ^( 2460 )“ 

1.31 ± 

0.04 ± 0.02 ± 

0.02 

1.31 ± 

0.04 ± 0.04 ± 

0.00 

24 ^( 2760 )“ 

0.053 + 

0.011 1 r )(-|0 1 

0.006 U.UUO =C 

0.008 

0.075 1 

aoo 8 ± 0.005 ± 

0.003 


8.4 Branching fractions 

The measured branching fraction of the —)■ decay in the phase-space region 

m(Z4°7r^) > 2.1 GeV/c^ is 

Z4°7r+7r“) = (8.46 ± 0.14 ± 0.29 ± 0.40) x 10“^ (33) 


taking into account the systematic uncertainties reported in Table The first uncertainty 
is statistical, the second systematic, and the third the uncertainty from the branching 
fraction of the B^ —)■ Z4*(2010)“7r’'' normalisation decay channel. The result agrees with 
the previous Belle measurement (8.4 ± 0.4 ± 0.8) x 10“^ 21 and the BaBar measurement 
(8.81 ± 0.18 ± 0.76 ± 0.78 ± 0.11)x 10“^ |^, obtained in a slightly larger phase-space 
region. A multiplicative factor of 94.5% (96.2%) is required to scale the Belle (BaBar) 
results to the same phase-space region as in this analysis. 
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Table 8: The phase of the complex coefficients of the resonant contributions for the Isobar 
model and the K-matrix model. The first uncertainty is statistical, the second and the third are 
experimental and model-dependent systematic uncertainties, respectively. 


Resonance 

Isobar (ari 

?(Ci)°) 


K-matrix (arg(cj)°' 

) 

Nonresonance 

77.1 ± 

4.5 T 

2.3 T 

5.4 


n/a 


/o(500) 

38.4 ± 

2.7 T 

1.3 T 

3.7 


n/a 


/o(980) 

138.9 ± 

4.6 T 

1.5 T 10.9 


n/a 


/o(2020) 

258.5 ± 

5.0 T 

1.1 T 26.8 


n/a 


p(770) 


0.0 (fixed) 



0.0 (fixed) 


w(782) 

176.8 ± 

7.8 T 

0.6 T 

0.5 

174.8 T 

8.0 T 1.5 T 

0.5 

p(1450) 

149.0 ± 

7.5 T 

4.8 T 

4.5 

132.9 T 

7.8 T 8.5 T 

5.5 

p(1700) 

103.5T13.it 

4.5 T 

2.4 

77.6 T 

9.9 T 23.1 T 

4.5 

/2(1270) 

158.1 T 

3.0 T 

1.6 T 

3.8 

147.8 T 

2.5 T 8.5 T 

2.6 

D^tt~ P-wave 

266.7 T 

3.7 T 

0.3 T 

7.1 

261.0 T 

4.0 T 3.3 T 

6.7 

T)(5(2400)- 

83.6 T 

4.4 T 

2.8 T 

4.6 

78.4 T 

4.1 T 11.5 T 

1.7 

Di{2Am)- 

262.9 T 

2.9 T 

0.8 T 

3.0 

257.4 T 

3.4 T 0.7 T 

1.9 

T)|(2760)- 

91.1 T 

6.7 T 

1.4 T 

5.1 

92.7 T 

7.3 T 15.2 T 

2.3 


Table 9: The fit fractions of the resonant contributions for the Isobar and K-matrix models 
with m{D^7r^) > 2.1 GeV/c^. The first uncertainty is statistical, the second and the third are 
experimental and model-dependent systematic uncertainties, respectively. 


Resonance Isobar %) K-matrix %) 


Nonresonance 

2.82 

T 

0.34 

T 

0.07 

T 

0.80 



n/a 



/o(500) 

13.2 

T 

0.89 

T 

0.31 

T 

2.45 



n/a 



/o(980) 

1.56 

T 

0.29 

T 

0.11 

T 

0.54 



n/a 



/o(2020) 

1.58 

T 

0.36 

T 

0.15 

T 

1.00 



n/a 



S-wave 

16.39 

T 

0.58 

T 

0.43 

T 

1.46 

16.51 

T 

0.70 T 1.68 

T 

1.10 

p(770) 

37.54 

T 

1.00 

T 

0.61 

T 

0.98 

36.15 

T 

1.00 T 2.13 

T 

0.79 

w(782) 

0.49 

T 

0.13 

T 

0.01 

T 

0.03 

0.50 

T 

0.13 T 0.01 

T 

0.02 

p(1450) 

1.54 

T 

0.32 

T 

0.08 

T 

0.22 

2.16 

T 

0.42 T 0.82 

T 

0.21 

p(1700) 

0.38 

-1- 

0.25 

0.12 

T 

0.07 

T 

0.06 

0.83 

T 

0.21 T 0.61 

T 

0.12 

/2(1270) 

10.28 

T 

0.49 

T 

0.31 

T 

1.10 

9.88 

T 

0.58 T 0.83 

T 

0.58 

D^tt~ P-wave 

9.21 

T 

0.56 

T 

0.24 

T 

1.73 

9.22 

T 

0.58 T 0.67 

T 

0.75 

L>S(2400)- 

9.00 

T 

0.60 

T 

0.20 

T 

0.35 

9.27 

T 

0.60 T 0.86 

T 

0.52 

L>^(2460)- 

28.83 

T 

0.69 

T 

0.74 

T 

0.50 

28.13 

T 

0.72 T 1.06 

T 

0.54 

L»^(2760)- 

1.22 

T 

0.19 

T 

0.07 

T 

0.09 

1.58 

T 

0.22 T 0.18 

T 

0.07 


The branching fraction of each qnasi-two-body decay, —)■ Vih^, with r* —)• hih 2 , is 

given by 

^ r,h) X B{r, ^ hihs) = B{B^ ^ D^7i+n-) x (34) 

where the resonant states {hih 2 ) = {D^7r~), (tt+tt"). The fit fractions T), defined in Eq. (|^, 
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Table 10: Correction factors due to the T>*(2010) veto. 


Resonance 

% 

/o(500) 

99.52 ±0.10 

/o(980) 

98.74 ±0.09 

/o(2020) 

99.29 ±0.05 

S-wave 

98.55 ±0.04 

p(770) 

98.95 ±0.03 

uj{782) 

99.39 ±0.02 

p(1450) 

95.66 ±0.06 

p(1700) 

96.73 ±0.06 

/2(1270) 

91.91 ±0.09 

Zl*(2400)- 

98.60 ±0.10 

Zl*(2460)- 

100. 

Zl*(2760)- 

100. 


are obtained from the Dalitz plot analysis and are listed in Table The correction factors, 
account for the cut-off due to the T)*(2010)“ veto. They are obtained by generating 
pseudo experiment samples for each resonance over the Dalitz plot and applying the same 
requirement (m(D°7r^) > 2.1 GeV/c^). They are summarised in Table 10 The correction 
factors are the same for the Isobar model and the K-matrix model. The effects due to the 
uncertainties of the masses and widths of the resonances are included in the uncertainties 
given in the table. 

Using the overall —)■ decay branching fraction, the fit fractions (J^j) 

and the correction factors the branching fractions of quasi two-body decays are 

calculated in Table [ll| The hrst observation of the decays B 0 ^ DO/o(500), B^ ^ 
D°/o(980), B^ —D°p(1450), as well as B^ —Zl3(2760)“7r+, and the hrst evidence 
of B^ — 11 *^/ 0 ( 2020 ) are reported. The present world averages 32 of the branching 
fractions B{B^ Wp{77Q)) x I3(p(770) ^ tt+tt-), B{B^ ^0/2(1270)) x i3(/2(1270) ^ 
TT+TT-), B{B^ Zl*(2400)-7r+) x i3(Zl*(2400)- ^ DOvt”), and B{B^ Zl*(2460)-7r+) x 
i3(Il*(2460)- ^ W-K-) are improved considerably. When accounting for the branching 
fractions of the a;(782) and /2(1270) to tt+tt”, one obtains the following results for the 
Isobar model 


B{B^ D°uj{782)) = (2.75 ± 0.72 ± 0.13 ± 0.20 ± 0.13/)J|5j) x 10”^ (35) 

and 

B{B^ DV2(1270)) = (16.8 ± 1.1 ± 0.7 ± 1.8 ± 0.7/|(;^) x 10”^ (36) 

For the K-matrix model, one obtains 

B{B^ Zi°a;(782)) = (2.81 ± 0.72 ± 0.13 ± 0.13 ± 0.13/)J;^^) x 10”^ (37) 

and 

B{B^ DV2(1270)) = (16.1 ± 1.1 ± 1.4 ± 0.9 ± 0.7/|(i) x 10”^ (38) 
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Table 11: Measured branching fractions of —)• rh^) x B{r —?■ /11/12) for the Isobar and 

K-matrix models. The first uncertainty is statistical, the second the experimental systematic, the 
third the model-dependent systematic, and the fourth the uncertainty from the normalisation 
—)■ L)*(2010)“7r''“ channel. 


Resonance Isobar (xlO K-matrix (xlO 


/o(500) 

11.2 ± 

0.8 ± 

0.5 ± 

2.1 ± 

0.5 




n/a 




fo(980) 

1.34 ± 

0.25 ± 

0.10 ± 

0.46 ± 

0.06 




n/a 




fo(2020) 

1.35 ± 

0.31 ± 

0.14 ± 

0.85 ± 

0.06 




n/a 




S-wave 

14.1 ± 

0.5 ± 

0.6 ± 

1.3 ± 

0.7 

14.2 ± 

0.6 

± 

1.5 

± 

0.9 ± 

0.7 

p(770) 

32.1 ± 

LOT 

1.2 ± 

0.9 ± 

1.5 

31.0 ± 

1.0 

± 

2.1 

± 

0.7 ± 

1.5 

a;(782) 

0.42 ± 

0.11 ± 

0.02 ± 

0.03 ± 

0.02 

-H 

CO 

o 

0.11 

± 

0.02 

± 

0.02 ± 

0.02 

p(1450) 

1.36 ± 

0.28 ± 

O 

O 

00 

H- 

0.19 ± 

0.06 

1.91 ± 

0.37 

± 

0.73 

± 

0.19 ± 

0.09 

p(1700) 

-H 

CO 

CO 

o 

0.11 ± 

0.06 ± 

0.05 ± 

0.02 

0.73 ± 

0.18 

± 

0.53 

± 

0.10 ± 

0.03 

/2(1270) 

9.5 ± 

0.5 ± 

0.4 ± 

1.0 ± 

0.4 

9.1 ± 

0.6 

± 

0.8 

± 

0.5 ± 

0.4 

L>(5(2400)- 

7.7 ± 

0.5 ± 

0.3 ± 

0.3 ± 

0.4 

8.0 ± 

0.5 

± 

0.8 

± 

0.4 ± 

0.4 

L>^(2460)- 

24.4 ± 

0.7 ± 

1.0 ± 

0.4 ± 

1.2 

23.8 ± 

0.7 

± 

1.2 

± 

0.5 ± 

1.1 

71^(2760)- 

1.03 ± 

0.16 ± 

-H 

o 

o 

-H 

00 

o 

o 

0.05 

1.34 ± 

0.19 

± 

0.16 

± 

0.06 ± 

0.06 


In both models, the hfth uncertainty is due to knowledge of the tt+tt decay rates j^. 
The results are consistent with the measurement of the decay —)■ D^u{782), using the 


dominant a;(782) —)■ tt+tt 7r° decay 32,40 


8.5 Structure of the /o(980) and /o(500) resonances 

In the Isobar model, signihcant contributions from both Z)°/o(500) and B^ —)• 

D^fo{980) decays are observed. The related branching fraction measurements can be used 
to obtain information on the substructure of the /o(980) and /o(500) resonances within the 
factorisation approximation. As discussed in Sec. two models for the quark structure of 
those states are considered: qq or [qq'][qq'] (tetraquarks). In both models, mixing angles 
between different quark states are determined using our measurements. In the qq model, 
the mixing between ss and uu or dd can be written as 

|/o(980)) = cos(^mix|ss)-fsin^Jmixl^h), (39) 

|/o(500)) = -sin<y9mix|ss)cos</9mix|nfi), (40) 

where \nn) = {\uu) + \dd))/\/2 and (^mix is the mixing angle. In the [qq'][qq'] model, the 
mixing angle, Wmix, is introduced and the mixing becomes 


|/o(980)) = cos WmixIrmss) -|- sincUmixIMwdd), 
|/o(500)) = — smu^i:i^\nnss) + cosojrnm\uudd). 

In both cases, the following variable is defined 


_ 


13{B^ T)°/o(980)) <I)(500) 
B{B0 BOfo(500)) ^ d>(980)’ 


(41) 

(42) 


(43) 
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Table 12: Systematic uncertainties on . The sum in quadrature of the uncertainties is also 
reported. 


Source 


BID 

0.001 

Trigger 

0.001 

Reconstruction 

< 0.001 

Simulation statistic 

0.001 

Background model 

0.001 

T>*(2010)- veto 

0.012 

Other res. 

0.007 

RBW parameters 

0.008 

TTTT res. mass, width 

0.011 

/o(500) model 

0.033 

/o(980) model 

0.028 

Total 

0.048 


where $(500) and $(980) are the integrals of the phase-space factors computed over 
the resonant lineshapes and the phase-space factors are proportional to the momentum 
computed in the rest frame. The value of their ratio is $(500)/$(980) = 1.02 ± 0.05. 

The value of the branching fraction i3(/o(500) — )■ 7r’''7r“) = 2/3 is obtained from the 
isospin Clebsch-Gordan coefficients and assumes that there are only contributions from tttt 
hnal states. The ratio i3(/o(980) —)■ K~)/B{fo{980) —)■ tt+tt") = 0 .35loj4, obtained 
from an average of the measurements by the BaBar and BES collaborations, is 
used to estimate the branching fraction B{fo{980) —)■ tt^tt”). Assuming that the tttt and 
KK decays are dominant in the /o(980) decays, B{fo{980) —)• tt+tt") = 0.46 ± 0.06 is 
obtained. This gives 

— a -I 77+0.066 

' — U.i/ <_o, 062 ) 


taking into account the systematic uncertainties, as listed in Table 12 
The parameter is related to the mixing angle by the equation 


= tan^ (pmix x 


F{B^ ^ /o(980)) 
F(50 ^ /o(500)) 


(44) 


in the qq model and by 


1 — v^tancjmix 

2 

F(R0 ^ /o(980)) 

tan cumix + \/2 


F(R0 ^ /o(500)) 


(45) 


in the tetraquark model 57,58 . The form factors F{B^ —)■ /o(980)) and F{B^ 


-)■ 


/o(500)) are evaluated at the four-momentum transfer squared equal to the square of 
the mass. Finally, values of the mixing angles as a function of form factor ratio are 
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Figure 13: Mixing angle as a function of form factor ratio for the (a) qq model and (b) [(/gild'll 
tetraquark model. Green band gives Icr interval around central values (black solid line). 
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49 


obtained in Fig. 13 for the qq model and the tetraquark model. Such angles have 

also been computed by LHCb for the decays —)■ J/ipn^'K~ 

The expectation is that the ratio of form factors should be 
LHCb has recently performed a search for the decay —)■ D^f q{980) 


chose 


to unity. However, 
The limit set 


95 


on this decay is below the value expected in a simple model based on our measured value 
of —)■ i4°/o(500)) and assuming equal form factors. More complicated models may 

be needed in order to explain all results. 

The above discussion is one possible interpretation of the results. Another possi¬ 
ble mechanism 91,96 involves the generation of pseudo-scalar resonances through the 
interactions of 7r'''7r“ mesons. 


11 


can 


8.6 Isospin analysis of the B Dp system 

The measured branching fraction of the B^ —)■ i7°p(770)° decay, presented in Table 
be used to perform an isospin analysis of the B ^ Dp system. Isospin symmetry relates 
the amplitudes of the decays —)■ T)°p(770)+, 5° —)■ D~p{770)^, and 5° —)■ T)°p(770)°, 

which can be written as linear combinations of the isospin eigenstates Aj with 1 = 1/2 
and 3/2 l^liT] 


leading to 


A{D°p+) 

— A/ 3 A 3 / 2 , 

(46) 

A{D-p+) 

= ^ 173 ^ 3/2 + ^2/3^1/2 5 


A{D^p^) 

= y^A3/2 - ^I73Ai/2, 


A{D^p+) 

= A{D-p+) + V2A{D^p^). 

(47) 


The strong phase difference between the amplitudes A 1/2 and A 3/2 is denoted by S^p- 
Final-state interactions between the states D^p^ and D~p^ may lead to a value of 5Dp 
different from zero and through constructive interference, to a larger value of B{B^ —)■ D^p^) 
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Table 13: 

Results of Rop 

and cos 5Dp- 

Model 

Rdp 

cos Sdp 

Isobar 

0.69 ±0.15 

u.aou_g Q4g 

K-matrix 

0.69 ±0.15 

u.ao (_o,o48 


than the prediction obtained within the factorisation approximation. In the heavy-quark 


limit, the factorisation model predicts 97,98 Sop = 0 {Aqcd/ fnb) and the amplitude ratio 
Rop = I = 1 -|- 0{AQCD/'mb), where represents the b quark mass and Aqcd the 


QCD scale. 

Using our measurement of D^p^) together with the world average values of 

13(5° —> 13(5+ —5°p+), and the ratio of lifetimes r(5+)/r(5°) 32 , we obtain 


Rdp — \l - 


3(5(5-p+) + 5(5V°)) 

- X - 


5(50p^ 


TbO 


1/2 


and 


cos 6dp = 


45 


X 


Dp 


Tb+ 

TbO 


X 


3(5(5-p+) -2 5(5V°)) 
5(5V) 


+ 1 


(48) 


(49) 


With a frequentist statistical approach 99 , Rop and cosbop are calculated for the Isobar 

These results are not signihcantly different from the 

and 


and K-matrix models in Table 13 


predictions of factorisation models. As opposed to the theoretical expectations 
in contrast to the 5 —system 


40 


37,41 


non-factorisable hnal-state interaction effects do 
not introduce a sizeable phase difference between the isospin amplitudes in the B ^ Dp 
system . The precision on Rop and cosSop is dominated by that of the branching fractions 


of the decays 5+ —5°p(770)+ (14%) and 5° —)■ 5 p(770)+ (17%) . The precision 

of the branching fraction of the 5° —)■ 5°p(770)° decay is 7.3% (9.2%) for the Isobar 


(K-matrix) model (see Table 11). 


9 Conclusion 

A Dalitz plot analysis of the 5° — 5°7r+7r“ decay is presented. The decay model contains 
four components from resonances, four P-wave tt+tt" resonances and one D-wave 

TT+TT” resonance. Two models are used to describe the S-wave tt+tt” resonances. The 
Isobar model uses four components, including the /o(500), /o(980), /o(2020) resonances 
and a nonresonant contribution. The K-matrix approach describes the tt+tt" S-wave using 
a 5 X 5 scattering matrix with a production vector. The overall branching fraction of 
5° —)■ 5°7r+7r“ and quasi-two-body decays are measured. Signihcant contributions from 
the /o(500), /o(980), p(1450) and 53(2760)“ mesons are observed for the hrst time. For 
the latter, this is a conhrmation of the observation from previous inclusive measurements, 
and the spin-parity of this resonance is determined for the hrst time to be = 3“. 
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This suggests a spectroscopic assignment of and shows that the ID family of charm 
resonances can be explored in Dalitz plot analysis of i?-meson decays in the same way as 
recently seen for the charm-strange resonances 25,26 . Evidence for the /o(2020) meson 


is also seen for the first time. The measured branching fractions of two-body decays are 
more precise than the existing world averages and there is good agreement between values 
from the Isobar and K-matrix models. 

The masses and widths of the resonances are also determined. The measnred 

masses and widths of the Zi)Q(2400)“ and ^3(2760)“ states are consistent with the previous 
measurements. The precision on the Dq(2400)~ meson is mnch improved. For the 
measnrement on the mass and width of the 172(2460)“ meson, the broad i7g(2400)“ 
component was excluded from the fit model in the former LHCb inclusive analysis due 
to a high correlation with the background lineshape parameters, while here it is included. 
The present result therefore supersedes the former measurement. 

The signihcant contribntions fonnd for both the /o(500) and /o(980) allow ns to 
constrain on the mixing angle between the /o(500) and /o(980) resonances. An isospin 
analysis in the B ^ Dp decays using our improved measurement of the branching fraction 
of the decay —)• D^p^ is performed, indicating that non-factorisable effects from 

hnal-state interactions are limited in the Dp system. 
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Appendices 


A Unnormalised Legendre polynomial weighted mo¬ 
ments 

Figures [14] and [T^ show the distributions of the unnormalised Legendre polynomial weighted 
moments < > which display the contributions of resonances with spin larger than 

L/2. The p(770) resonance can clearly be seen in the distributions with L <2 and the 
i72(2460)“ resonance in the distributions with L < 4. Figures 16 and 17 display an 


expanded version in low mass regions. The distributions from the Isobar and the K-matrix 
models are compatible with those from data. 
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Figure 14: The first eight unnormalised Legendre polynomial weighted moments (0 to 7 correspond 
to (a) to (h)) for background-subtracted and efficiency-corrected > I1^7r'*'7r“ data and the 
Dalitz plot fit results as a function of m'^(D^Tr~). 
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Figure 15: The first eight unnormalised Legendre polynomial weighted moments (0 to 7 correspond 
to (a) to (h)) for background-subtracted and efficiency-corrected —?■ I1^7r'*'7r“ data and the 
Dalitz plot fit results as a function of m?' k~). 
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Figure 16: The first eight unnormalised Legendre polynomial weighted moments (0 to 7 correspond 
to (a) to (h)) for background-subtracted and efficiency-corrected —?■ data and the 

Dalitz plot fit results as a function of Only results in the region < 10 

GeV^/c^ are shown. 
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Figure 17: The first eight unnormalised Legendre polynomial weighted moments (0 to 7 correspond 
to (a) to (h)) for background-subtracted and efficiency-corrected —?■ I1^7r'’'7r“ data and the 
Dalitz plot fit results as a function of m^(7r'*“7r“). Only results in the region m‘^(7r~^7r~) < 3 
GeV^/c^ are shown. 
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B Systematic uncertainties on the parameters in the 
Dalitz plot analysis 

B.l Systematic uncertainties for the Isobar model 


Table 14: Systematic uncertainties on the D^'k resonant masses (MeV/c^) and widths (MeV) 
for the Isobar model. 


Source 

£>5(2400)- 

£5(2460)- 

£5(2760)- 

To 

mo 

To 

mo 

To 

mo 

PID 

1.9 

0.5 

< 0.1 

< 0.1 

1.2 

0.3 

Trigger 

0.5 

0.2 

0.2 

< 0.1 

1.9 

0.7 

Reconstruction 

0.2 

0.1 

< 0.1 

< 0.1 

0.5 

0.1 

Simulation statistic 

0.6 

0.1 

0.1 

< 0.1 

0.6 

0.1 

Background model 

1.5 

0.6 

0.1 

< 0.1 

3.4 

0.4 

£>*(2010)- veto 

4.4 

0.8 

< 0.1 

< 0.1 

4.6 

0.4 

Total (experiment) 

5.1 

1.1 

0.3 

< 0.1 

6.2 

1.0 

Additional resonances 

10.7 

0.4 

0.1 

0.1 

21.0 

5.1 

RBW parameters 

0.1 

1.9 

0.1 

< 0.1 

1.0 

1.5 

TTTT res. mass, width 

3.4 

0.1 

0.1 

0.1 

2.9 

1.3 

B^, mass 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

/o(500) model 

2.3 

3.6 

0.5 

< 0.1 

9.0 

3.5 

/o(980) model 

2.8 

0.7 

0.2 

0.1 

3.5 

1.2 

Total (model) 

11.8 

4.2 

0.6 

0.3 

23.3 

6.6 

Total (all) 

12.9 

4.3 

0.7 

0.3 

24.1 

6.7 
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Table 15: Systematic uncertainties on the moduli of the complex coefficients of the resonant 
contributions for the Isobar model. The moduli are normalised to that of p{770). 


Source 

Nonres. 

/o(500) 

a;(782) 

/o(980) 

/2(1270) 

p(1450) 

PID 

0.02 

0.15 

< 0.01 

0.03 

< O.OOl 

0.01 

Trigger 

0.02 

0.03 

< 0.01 

0.01 

< 0.001 

< 0.01 

Reconstruction 

< 0.01 

0.01 

< 0.01 

< 0.01 

< 0.001 

< 0.01 

Simulation statistic 

< 0.01 

0.11 

< 0.01 

0.01 

< 0.001 

< 0.01 

Background model 

< 0.01 

0.07 

< 0.01 

0.02 

< 0.001 

< 0.01 

D*(2010)- veto 

0.03 

0.20 

< 0.01 

0.08 

0.002 

0.01 

Total (experiment) 

0.04 

0.29 

< 0.01 

0.09 

< 0.001 

0.01 

Additional resonances 

0.34 

0.61 

< 0.01 

0.03 

0.003 

< 0.01 

RBW parameters 

0.11 

0.10 

0.01 

0.04 

0.001 

0.01 

TTTT res. mass, width 

0.06 

0.41 

< 0.01 

0.03 

0.001 

< 0.01 

/o(500) model 

0.36 

n/a 

< 0.01 

0.46 

0.004 

0.01 

/o(980) model 

0.01 

0.30 

< 0.01 

n/a 

< 0.001 

0.01 

Total (model) 

0.51 

0.80 

0.01 

0.46 

0.005 

0.02 

Total 

0.51 

0.85 

0.01 

0.47 

0.006 

0.02 


Source 

p(1700) 

/o(2020) 

P-wave 

£>5(2400)- 

£>5(2460)- 

£>1(2760)- 

PID 

< 0.001 

0.06 

0.14 

0.21 

< 0.01 

< 0.001 

Trigger 

0.001 

0.02 

0.17 

0.08 

0.02 

< 0.001 

Reconstruction 

< 0.001 

< 0.01 

0.01 

0.02 

< 0.01 

< 0.001 

Simulation statistic 

0.001 

0.02 

0.10 

0.03 

0.01 

< 0.001 

Background model 

0.002 

0.01 

0.03 

0.06 

< 0.01 

0.002 

i:)*(2010)- veto 

0.006 

0.20 

n/a 

0.20 

< 0.01 

0.002 

Total (experiment) 

0.006 

0.21 

0.25 

0.31 

0.02 

0.003 

Additional resonances 

0.005 

0.37 

1.74 

0.47 

0.01 

0.007 

RBW parameters 

0.002 

0.27 

0.05 

0.12 

0.01 

< 0.001 

TTTT res. mass, width 

0.002 

0.73 

0.16 

0.19 

< 0.01 

0.001 

/o(500) model 

0.004 

1.56 

0.62 

0.20 

0.01 

0.003 

/o(980) model 

0.003 

0.06 

0.06 

0.18 

0.01 

0.002 

Total (model) 

0.008 

1.78 

1.86 

0.59 

0.02 

0.008 

Total (all) 

0.010 

1.80 

1.87 

0.66 

0.03 

0.008 
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Table 16; Systematic uncertainties on the phases (°) of the complex coefficients of the resonant 
contributions for the Isobar model. The phase of /9(700) is set to 0° as the reference. 


Source 

Nonres. 

/o(500) 

a;(782) 

/o(980) 

/2(1270) 

p(1450) 

PID 

0.2 

0.1 

0.4 

0.1 

0.1 

0.1 

Trigger 

1.1 

0.3 

0.1 

0.1 

0.6 

0.2 

Reconstruction 

< 0.1 

< 0.1 

< 0.1 

< 0.1 

0.1 

0.1 

Simulation statistic 

0.3 

0.1 

0.3 

0.2 

0.3 

0.3 

Background model 

0.7 

0.1 

< 0.1 

0.4 

0.4 

0.4 

£>*(2010)- veto 

1.8 

1.2 

0.2 

1.4 

1.4 

4.8 

Total (experiment) 

2.3 

1.3 

0.6 

1.5 

1.6 

4.8 

Additional resonances 

3.0 

2.8 

< 0.1 

2.7 

1.6 

1.1 

RBW parameters 

0.3 

1.1 

0.1 

1.3 

< 0.1 

0.3 

TTTT res. mass, width 

2.3 

0.1 

0.3 

0.4 

1.2 

2.2 

/o(500) model 

3.7 

n/a 

0.3 

10.5 

3.0 

3.6 

/o(980) model 

1.0 

2.1 

0.2 

n/a 

1.1 

1.1 

Total (model) 

5.4 

3.7 

0.5 

10.9 

3.8 

4.5 

Total (all) 

5.8 

3.9 

0.7 

11.0 

4.1 

6.6 


Source 

p(1700) 

/o(2020) 

D^tt P-wave 

£>5(2400)- 

£5(2460)- 

£|(2760)- 

PID 

0.4 

< 0.1 

0.2 

0.3 

0.1 

0.2 

Trigger 

0.5 

0.2 

0.2 

0.7 

< 0.1 

0.6 

Reconstruction 

< 0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

Simulation statistic 

0.7 

0.2 

0.1 

0.2 

0.1 

0.1 

Background model 

0.6 

0.9 

0.1 

< 0.1 

0.2 

0.3 

£>*(2010)- veto 

4.4 

0.6 

n/a 

2.7 

0.8 

1.2 

Total (experiment) 

4.5 

1.1 

0.3 

2.8 

0.8 

1.4 

Additional resonances 

0.3 

3.4 

5.7 

1.8 

1.2 

3.9 

RBW parameters 

1.0 

1.1 

0.6 

1.6 

0.7 

1.8 

TTTT res. mass, width 

0.5 

1.3 

0.5 

0.3 

0.7 

2.0 

/o(500) model 

0.3 

26.5 

4.2 

3.9 

2.4 

1.9 

/o(980) model 

2.1 

1.8 

0.4 

0.1 

0.9 

0.2 

Total (model) 

2.4 

26.8 

7.1 

4.6 

3.0 

5.1 

Total (all) 

5.1 

26.9 

7.1 

5.4 

3.1 

5.3 
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Table 17: Systematic uncertainties on the fit fractions (%) of the resonant contributions for the 
Isobar model. 


Source 

Nonres. 

/o(500) 

p(770) 

a;(782) 

/o(980) 

/2(1270) 

p(1450) 

PID 

0.02 

0.23 

0.37 

0.01 

0.04 

0.07 

0.06 

Trigger 

0.01 

0.03 

0.09 

< 0.01 

< 0.01 

0.05 

0.03 

Reconstruction 

0.01 

0.09 

0.39 

0.01 

0.01 

0.01 

< 0.01 

Simulation statistic 

< 0.01 

0.09 

0.26 

< 0.01 

0.01 

0.04 

< 0.01 

Background model 

0.01 

0.08 

0.07 

< 0.01 

0.01 

0.05 

0.04 

_D*(2010)- veto 

0.06 

0.15 

0.01 

< 0.01 

0.10 

0.29 

0.03 

Total (experiment) 

0.07 

0.31 

0.61 

0.01 

0.11 

0.31 

0.08 

Additional resonances 

0.51 

0.57 

0.81 

0.02 

< 0.01 

0.53 

0.01 

RBW parameters 

0.21 

0.03 

0.50 

0.01 

0.07 

0.10 

0.15 

TTTT res. mass, width 

0.09 

0.60 

0.08 

< 0.01 

0.03 

0.16 

< 0.01 

/o(500) model 

0.57 

2.25 

0.04 

0.01 

0.50 

0.94 

0.14 

/o(980) model 

0.03 

0.48 

0.19 

0.01 

0.20 

0.11 

0.09 

Total (model) 

0.80 

2.45 

0.98 

0.03 

0.54 

1.10 

0.22 

Total (all) 

0.80 

2.47 

1.15 

0.03 

0.56 

1.14 

0.24 


Source 

p(1700) 

/o(2020) 

S-wave 

P-wave 

D*(2400)- 

D5(2460)- 

D*(2760)- 

PID 

0.01 

0.01 

0.31 

0.19 

0.12 

0.43 

0.02 

Trigger 

< 0.01 

0.03 

0.04 

0.02 

0.06 

0.28 

0.01 

Reconstruction 

< 0.01 

0.01 

0.12 

0.13 

0.14 

0.18 

0.02 

Simulation statistic 

< 0.01 

0.01 

0.05 

0.06 

0.04 

0.09 

0.01 

Background model 

0.02 

0.01 

0.04 

0.05 

0.03 

0.01 

0.02 

D*(2010)- veto 

0.07 

0.14 

0.26 

n/a 

0.01 

0.50 

0.06 

Total (experiment) 

0.07 

0.15 

0.43 

0.24 

0.20 

0.74 

0.07 

Additional resonances 

0.04 

0.24 

1.43 

1.60 

0.01 

0.01 

0.07 

RBW parameters 

0.02 

0.21 

0.07 

0.17 

0.28 

0.44 

0.02 

TTTT res. mass, width 

0.02 

0.22 

0.05 

0.17 

0.13 

0.02 

0.02 

/o(500) model 

0.03 

0.92 

0.25 

0.60 

0.13 

0.22 

0.03 

/o(980) model 

0.02 

0.04 

0.14 

0.01 

0.10 

0.09 

0.03 

Total (model) 

0.06 

1.00 

1.46 

1.73 

0.35 

0.50 

0.09 

Total (all) 

0.10 

1.01 

1.52 

1.74 

0.40 

0.90 

0.11 
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B.2 Systematic uncertainties for the K-matrix model 


Table 18: Systematic uncertainties on the D^tt resonant masses (MeV/c^) and widths (MeV) 
for the K-matrix model. 


Source 

£>5(2400)- 

£>2(2460)- 

£>5(2760)- 

To 

Too 

To 

mo 

To 

mo 

PID 

8.9 

5.9 

1.0 

0.3 

1.1 

8.6 

Trigger 

1.1 

0.6 

0.1 

< 0.1 

0.1 

0.1 

Reconstruction 

< 0.1 

< 0.1 

< 0.1 

< 0.1 

< 0.1 

< 0.1 

Simulation statistic 

4.6 

0.5 

0.2 

0.1 

2.4 

0.4 

Background model 

2.3 

0.9 

0.1 

< 0.1 

2.7 

0.8 

D*{2010)- veto 

14.9 

9.0 

1.4 

0.3 

12.1 

5.4 

Total (experiment) 

18.1 

10.8 

1.7 

0.4 

12.7 

10.2 

Additional resonances 

8.5 

1.0 

0.3 

0.2 

1.4 

1.1 

RBW parameters 

6.0 

1.5 

< 0.1 

0.1 

8.9 

2.3 

TTTT res. mass, width 

1.4 

0.9 

0.1 

0.1 

1.3 

0.1 

mass 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

Total (model) 

10.5 

2.0 

0.4 

0.3 

9.1 

2.6 

Total (all) 

21.0 

11.0 

1.8 

0.5 

15.6 

10.5 
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Table 19: Systematic uncertainties on the moduli of the complex coefficients of the resonant 
contributions for the K-matrix model. The moduli are normalised to that of p{770). 


Source 

u;(782) 

/2(1270) 

/9(1450) 

p(1700) 

PID 

< 0.01 

0.001 

0.04 

0.037 

Trigger 

< 0.01 

< 0.001 

< 0.01 

0.002 

Reconstruction 

< 0.01 

< 0.001 

< 0.01 

< 0.001 

Simulation statistic 

< 0.01 

0.003 

0.02 

0.004 

Background model 

< 0.01 

< 0.001 

< 0.01 

0.001 

£'*(2010)- veto 

0.01 

0.005 

0.07 

0.067 

Total (experiment) 

0.01 

0.006 

0.08 

0.077 

Additional resonances 

< 0.01 

0.003 

< 0.01 

0.003 

RBW parameters 

0.01 

< 0.001 

0.01 

0.003 

TTTT res. mass, width 

< 0.01 

< 0.001 

< 0.01 

0.010 

Total (model) 

0.01 

0.003 

0.01 

0.011 

Total (all) 

0.01 

0.007 

0.08 

0.077 


Source 

U^TT P-wave 

D*(2400)- 

£>*(2460)- 

£>5(2760)- 

PID 

0.21 

0.10 

< 0.01 

0.003 

Trigger 

0.12 

0.01 

0.02 

< 0.001 

Reconstruction 

< 0.01 

< 0.01 

< 0.01 

0.001 

Simulation statistic 

0.62 

0.72 

0.03 

0.003 

Background model 

0.13 

0.13 

< 0.01 

0.002 

D*(2010)- veto 

< 0.01 

0.39 

< 0.01 

0.001 

Total (experiment) 

0.68 

0.84 

0.04 

0.005 

Additional resonances 

0.52 

0.22 

< 0.01 

0.002 

RBW parameters 

0.30 

0.46 

< 0.01 

0.002 

TTTT res. mass, width 

0.05 

0.04 

< 0.01 

0.001 

Total (model) 

0.60 

0.51 

< 0.01 

0.003 

Total (all) 

0.91 

0.98 

0.04 

0.006 
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Table 20; Systematic uncertainties on the phases (°) of the complex coefficients of the resonant 
contributions for the K-matrix model. The phase of /3(700) is set to 0° as reference. 


Source 

u;(782) 

/2(1270) 

/9(1450) 

p(1700) 

PID 

0.8 

5.6 

4.1 

12.8 

Trigger 

< 0.1 

0.2 

< 0.1 

0.3 

Reconstruction 

0.1 

< 0.1 

< 0.1 

0.0 

Simulation statistic 

0.5 

1.1 

1.6 

1.5 

Background model 

0.2 

1.3 

0.3 

0.1 

£•*(2010)- veto 

1.1 

6.2 

7.1 

19.2 

Total (experiment) 

1.5 

8.5 

8.4 

23.1 

Additional resonances 

0.2 

1.9 

2.3 

3.3 

RBW parameters 

0.4 

1.3 

0.4 

0.6 

TTTT res. mass, width 

0.1 

1.1 

5.0 

3.0 

Total (model) 

0.5 

2.6 

5.5 

4.5 

Total (all) 

1.5 

8.9 

10.0 

23.6 


Source 

D'V P-wave 

D*(2400)- 

£>*(2460)- 

£>5(2760)- 

PID 

3.2 

7.0 

0.2 

11.8 

Trigger 

0.4 

0.9 

0.3 

0.1 

Reconstruction 

< 0.1 

< 0.1 

< 0.1 

0.1 

Simulation statistic 

0.1 

0.8 

0.1 

1.3 

Background model 

0.7 

0.3 

0.4 

0.1 

D*(2010)- veto 

n/a 

9.1 

0.4 

9.5 

Total (experiment) 

3.3 

11.5 

0.7 

15.2 

Additional resonances 

6.6 

0.2 

1.8 

0.4 

RBW parameters 

1.0 

1.5 

0.4 

2.2 

TTTT res. mass, width 

0.4 

0.7 

0.2 

0.3 

Total (model) 

6.7 

1.7 

1.9 

2.3 

Total (all) 

7.5 

11.7 

2.0 

15.4 
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Table 21: Systematic uncertainties on the fit fractions (%) of the resonant contributions for the 
K-matrix model. 


Source 

p(770) 

w(782) 

/2(1270) 

p(1450) 

/9(1700) 

PID 

1.36 

< 0.01 

0.06 

0.50 

0.37 

Trigger 

0.33 

0.01 

0.04 

0.01 

0.01 

Reconstruction 

0.02 

< 0.01 

< 0.01 

< 0.01 

< 0.01 

Simulation statistic 

1.17 

0.01 

0.45 

0.22 

0.02 

Background model 

0.07 

< 0.01 

0.15 

0.04 

0.01 

£>*(2010)- veto 

1.09 

< 0.01 

0.68 

0.61 

0.48 

Total (experiment) 

2.13 

0.01 

0.83 

0.82 

0.61 

Additional resonances 

0.61 

0.02 

0.56 

0.01 

0.05 

RBW parameters 

0.49 

0.01 

0.15 

0.21 

0.04 

TTTT res. mass, width 

0.12 

0.01 

0.04 

0.04 

0.10 

Total (model) 

0.79 

0.02 

0.58 

0.21 

0.12 

Total (all) 

2.28 

0.03 

1.02 

0.85 

0.62 


Source 

S-wave 

P-wave 

£•5(2400)- 

£•5(2460)- 

£•5(2760)- 

PID 

0.77 

0.59 

0.13 

0.28 

0.13 

Trigger 

0.09 

0.03 

0.01 

0.39 

0.03 

Reconstruction 

0.01 

< 0.01 

< 0.01 

< 0.01 

< 0.01 

Simulation statistic 

0.31 

0.28 

0.49 

0.38 

0.01 

Background model 

0.23 

0.13 

0.08 

0.03 

0.04 

£>*(2010)- veto 

1.44 

n/a 

0.69 

0.86 

0.12 

Total (experiment) 

1.68 

0.67 

0.86 

1.06 

0.18 

Additional resonances 

1.08 

0.63 

0.23 

0.28 

0.03 

RBW parameters 

0.18 

0.40 

0.47 

0.46 

0.06 

TTTT res. mass, width 

0.07 

0.07 

0.02 

0.01 

0.01 

Total (model) 

1.10 

0.75 

0.52 

0.54 

0.07 

Total (all) 

2.01 

1.00 

1.01 

1.19 

0.20 
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C Results for the interference fit fractions 


The central values of the interference fit fractions for the Isobar (K-matrix) model are given 
in Table (Table [2^ . The statistical, experimental systematic and model-dependent 
uncertainties on these quantities are given in Tables 24 25 and 26 (Tables 27, 28 and 29). 


Table 22: Interference fit fractions (%) of the resonant contributions for the Isobar model with 
> 2.1 GeV/c^. The resonances are: (Tq) nonresonant S-wave, (Ti) /o(500), {A 2 ) 
M980), (T 3 ) /o(2020), (^ 4 ) p(770), (A) a;(782), [Aq) p(1450), (T 7 ) p(1700), (Ag) /2(1270), {Ag) 
D^7r~ P-wave, (^lo) i7g(2400)“, (^n) 772(2460)“, {A 12 ) 773(2760)“. The diagonal elements 
correspond to the fit fractions given in Table 



^0 

Ai 

Ag 

As 

A4 

As 

Aq 

Ar 

^8 

A9 

Aio 

All 

Ai2 

Ao 

2.82 

2.70 

-0.37 

-0.25 

0.00 

0.00 

0.00 

0.00 

0.13 

0.79 

-1.70 

-2.12 

0.06 

Ai 

- 

13.23 

-1.02 

-4.53 

0.00 

0.00 

0.00 

0.00 

0.14 

3.37 

0.97 

3.81 

0.57 

^2 

- 

- 

1.56 

0.79 

0.00 

0.00 

0.00 

0.00 

-0.21 

-0.60 

0.63 

-0.90 

-0.14 

As 

- 

- 

- 

1.58 

0.00 

0.00 

0.00 

0.00 

-0.17 

-1.39 

-1.27 

-1.76 

-0.16 

A4 

- 

- 

- 

- 

37.54 

-0.78 

2.43 

1.53 

0.00 

-5.71 

-1.54 

-3.26 

-0.78 

As 

- 

- 

- 

- 

- 

0.49 

-0.01 

0.00 

0.00 

0.00 

0.01 

-0.01 

0.00 

Aq 

- 

- 

- 

- 

- 

- 

1.54 

-0.06 

0.00 

0.26 

-0.74 

0.94 

0.04 

A7 

- 

- 

- 

- 

- 

- 

- 

0.38 

0.00 

-0.89 

-0.66 

-0.58 

-0.13 

As 

- 

- 

- 

- 

- 

- 

- 

- 

10.28 

-2.29 

-0.89 

-1.43 

-0.27 

Ag 

- 

- 

- 

- 

- 

- 

- 

- 

- 

9.21 

0.00 

-0.01 

0.00 

Aio 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

9.00 

0.01 

0.00 

All 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

28.83 

0.00 

Ai2 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

1.22 


Table 23: Interference fit fractions (%) of the resonant contributions for the K-matrix model with 
m( 77 ° 7 r=*") > 2.1 GeV/c^. The resonances are: ( 2 I 0 ) K-matrix S-wave, (Ai) p(770), (^ 2 ) a;(782), 
(T 3 ) p(1450), ( 7 I 4 ) p(1700), (^ 5 ) /2(1270), (^e) ^°vr“ P-wave, {Ay) 77*(2400)“, (As) 77^(2460)“, 
(Ag) 773 ( 2760 )“ The diagonal elements correspond to the fit fractions given in Table pll 



Aq 

Ai 

A2 

^3 

A4 

Aq 

Aq 

Ay 

^8 

^9 

^0 

16.51 

0.00 

0.00 

0.00 

0.00 

-0.06 

2.37 

-1.45 

-0.10 

0.01 


— 

36.15 

-0.84 

4.20 

2.10 

0.00 

-5.39 

-1.88 

-2.81 

-0.90 

A2 

— 

— 

0.50 

-0.01 

0.00 

0.00 

0.00 

0.01 

-0.01 

0.00 

A3 

— 

— 

— 

2.16 

-0.43 

0.00 

-0.15 

-1.14 

0.73 

-0.04 

A4 

— 

— 

— 

— 

0.83 

0.00 

-1.49 

-0.99 

-1.12 

-0.24 

A3 

— 

— 

— 

— 

— 

9.88 

-2.03 

-0.73 

-1.50 

-0.35 

Aq 

— 

— 

— 

— 

— 

— 

9.22 

0.00 

-0.01 

0.00 

Ay 

— 

— 

— 

— 

— 

— 

— 

9.27 

0.01 

0.00 

Ag 

- 

- 

- 

- 

- 

- 

- 

- 

28.13 

0.00 

Ag 

- 

- 

- 

- 

- 

- 

- 

- 

- 

1.58 
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Table 24: Statistical uncertainties on the interference fit fractions (%) of the resonant contribu¬ 
tions for the Isobar model with m(-D®7r^) > 2.1 GeV/c^. The resonances are: (^o) nonresonant 
S-wave, (Ai) /o(500), {A 2 ) /o(980)^(A3) /o(2020), (^ 4 ) p{770), (A) a;(782), (Ag) p(1450), 
(Aj) p(1700), (As) /2(1270), (Ag) D^n- P-wave, (Aio) D5(2400)-, (An) T»^(2460)-, (An) 
1 ) 3 ( 2760 )“. The diagonal elements correspond to the statistical uncertainties on the fit fractions 
given in Table 



Ao 

Ai 

A 2 

A 3 

A 4 

As 

Ag 

A 7 

As 

Ag 

Aio 

All 

Ai2 

^0 

0.34 

0.29 

0.11 

0.07 

0.00 

0.00 

0.00 

0.00 

0.02 

0.13 

0.36 

0.26 

0.03 

Ai 

— 

0.89 

0.54 

0.64 

0.00 

0.00 

0.00 

0.00 

0.04 

0.22 

0.45 

0.20 

0.08 

A 2 

— 

— 

0.29 

0.15 

0.00 

0.00 

0.00 

0.00 

0.03 

0.09 

0.09 

0.11 

0.04 

A 3 

— 

— 

— 

0.36 

0.00 

0.00 

0.00 

0.00 

0.02 

0.18 

0.20 

0.22 

0.05 

A 4 

— 

— 

— 

— 

1.00 

0.33 

0.65 

0.32 

0.00 

0.35 

0.21 

0.22 

0.12 


— 

— 

— 

— 

— 

0.13 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Ag 

— 

— 

— 

— 

— 

— 

0.32 

0.24 

0.00 

0.24 

0.18 

0.19 

0.04 

A 7 

— 

— 

— 

— 

— 

— 

— 

+0.25 

-0.12 

0.00 

0.23 

0.15 

0.20 

0.03 

^8 

— 

— 

— 

— 

— 

— 

— 


0.49 

0.20 

0.14 

0.07 

0.04 

Ag 

— 

— 

— 

— 

— 

— 

— 

— 

— 

0.56 

0.00 

0.00 

0.00 

Aio 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

0.60 

0.00 

0.00 

All 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0.69 

0.00 

Ai2 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

0.19 


Table 25: Experimental systematic uncertainties on the interference fit fractions (%) of the 
resonant contributions for the Isobar model with m{D^Tr^) > 2.1 GeV/c^. The resonances 
are: (Aq) nonresonant S-wave, (Ai) /o(500), (A 2 ) /o(980), (A 3 ) /o(2020), (A 4 ) p{770), (A 5 ) 
a;(782), (Ag) p(1450), (A 7 ) p(1700), (Ag) /2(1270), (Ag) D^n- P-wave, (Aio) 7)^(2400)-, (An) 
1 ) 2 ( 2460 )“, (A 12 ) 7 ) 3 ( 2760 )“. The diagonal elements correspond to the statistical uncertainties 
on the fit fractions given in Table 



Ao 

Ai 

A2 

A3 

A4 

Ag 

Ag 

A7 

As 

Ag 

Aio 

All 

Ai 2 

Ag 

0.07 

0.06 

0.03 

0.03 

0.00 

0.00 

0.00 

0.00 

0.07 

0.02 

0.11 

0.09 

0.01 

Ai 

— 

0.31 

0.22 

0.20 

0.00 

0.00 

0.00 

0.00 

0.07 

0.06 

0.26 

0.04 

0.01 

A2 

— 

— 

0.11 

0.05 

0.00 

0.00 

0.00 

0.00 

0.14 

0.02 

0.04 

0.04 

0.02 

A3 

— 

— 

— 

0.15 

0.00 

0.00 

0.00 

0.00 

0.11 

0.08 

0.06 

0.11 

0.02 

A4 

— 

— 

— 

— 

0.61 

0.05 

0.31 

0.16 

0.00 

0.18 

0.05 

0.03 

0.02 

Ag 

— 

— 

— 

— 

— 

0.01 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Ag 

— 

— 

— 

— 

— 

— 

0.08 

0.03 

0.00 

0.11 

0.12 

0.06 

0.01 

Ar 

— 

— 

— 

— 

— 

— 

— 

0.07 

0.00 

0.10 

0.07 

0.03 

0.01 

As 

— 

— 

— 

— 

— 

— 

— 

— 

0.31 

0.16 

0.11 

0.07 

0.02 

Ag 

— 

— 

— 

— 

— 

— 

— 

— 

— 

0.24 

0.00 

0.00 

0.00 

Aio 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

0.20 

0.00 

0.00 

All 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0.74 

0.00 

Ai2 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

0.07 
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Table 26: Model-dependent systematic uncertainties on the interference fit fractions (%) of 
the resonant contributions for the Isobar model with m{D^TT^) > 2.1 GeV/c^. The resonances 
are: (Aq) nonresonant S-wave, (Ai) /o(500), {A 2 ) fo{980), {A 3 ) /o(2020), (A 4 ) p{770), (A 5 ) 
a;(782), (Ag) p(1450), (A7) p(1700), (Ag) /2(1270), (Ag) D^tt- P-wave, (Aio) 77^(2400)-, (An) 
Z?2(2460)“, (A 12 ) 743 ( 2760 )“. The diagonal elements correspond to the statistical uncertainties 
on the fit fractions given in Table 



Aq 

Ai 

A 2 

A 3 

A 4 

As 

Aq 

Ay 

As 

Ag 

Aio 

All 

Ai2 

Ag 

0.80 

0.61 

0.31 

0.17 

0.00 

0.00 

0.00 

0.00 

0.03 

0.28 

0.56 

0.45 

0.01 

Ai 

— 

2.45 

2.00 

3.03 

0.00 

0.00 

0.00 

0.00 

0.16 

0.72 

0.79 

0.98 

0.08 

^2 

— 

— 

0.54 

0.67 

0.00 

0.00 

0.00 

0.00 

0.02 

0.15 

0.13 

0.28 

0.08 

A 3 

— 

— 

— 

1.00 

0.00 

0.00 

0.00 

0.00 

0.08 

0.51 

0.34 

0.69 

0.06 


— 

— 

— 

— 

0.98 

0.03 

0.47 

0.12 

0.00 

0.54 

0.33 

0.27 

0.09 

A 5 

— 

— 

— 

— 

— 

0.03 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Ag 

— 

— 

— 

— 

— 

— 

0.22 

0.08 

0.00 

0.31 

0.18 

0.14 

0.03 

Ay 

— 

— 

— 

— 

— 

— 

— 

0.06 

0.00 

0.12 

0.07 

0.04 

0.02 

As 

— 

— 

— 

— 

— 

— 

— 

— 

1.10 

0.49 

0.33 

0.09 

0.05 

Ag 

— 

— 

— 

— 

— 

— 

— 

— 

— 

1.73 

0.00 

0.00 

0.00 

Aio 











0.35 

0.00 

0.00 

All 












0.50 

0.00 

Ai2 













0.09 


Table 27: Statistical uncertainties on the interference fit fractions (%) of the resonant contribu¬ 
tions for the K-matrix model with m{D^TT^) > 2.1 GeV/c^. The resonances are: (Ao)K-matrix 
S-wave, (Ai) p(770), (A 2 ) a;(782), (A 3 ) p(1450), (A 4 ) /o(1700), (A 5 ) /2(1270), (Ag) 74 ° 7 r“ P-wave, 
(Ay) 74g(2400)“, (Ag) 742 ( 2460 )“, (Ag) 743(2760)“ The diagonal elements correspond to the 
statistical uncertainties on the fit fractions shown in Table 



Aq 

Ai 

A 2 

A 3 

A 4 

As 

Ag 

Ay 

As 

Ag 

^0 

0.70 

0.00 

0.00 

0.00 

0.00 

0.04 

0.28 

0.49 

0.43 

0.16 

^1 

— 

1.00 

0.34 

0.71 

0.31 

0.00 

0.41 

0.22 

0.25 

0.14 

A 2 

— 

— 

0.13 

0.00 

0.00 

0.00 

0.01 

0.00 

0.00 

0.00 

A 3 

— 

— 

— 

0.42 

0.29 

0.00 

0.29 

0.23 

0.20 

0.07 

A 4 

— 

— 

— 

— 

0.21 

0.00 

0.23 

0.16 

0.19 

0.04 

^5 

— 

— 

— 

— 

— 

0.58 

0.22 

0.16 

0.08 

0.05 

^6 

— 

— 

— 

— 

— 

— 

0.58 

0.00 

0.00 

0.00 

Ay 

— 

— 

— 

— 

— 

— 

— 

0.60 

0.00 

0.00 

As 

— 

— 

— 

— 

— 

— 

— 

— 

0.72 

0.00 

Ag 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0.22 
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Table 28 : Experimental systematic uncertainties on the interference fit fractions (%) of the 
resonant contributions for the K-matrix model with > 2.1 GeV/c^. The resonances 

are: (^o) K-matrix S-wave, (^i) / 9 ( 770 ), (^2) w( 782 ), (^43) / 9 ( 1450 ), (A 4 ) / 3 ( 1700 ), (^ 5 ) / 2 ( 1270 ), 
{Aq) D^'k~ P-wave, (^7) T)g( 2400 )“, {A^) 742 ( 2460 )“, (^9) 743 ( 2760 )“ The diagonal elements 
correspond to the statistical uncertainties on the fit fractions shown in Table 



Ao 

Ai 

A 2 

A 3 

A 4 

As 

As 

A 7 

As 

A^ 

^0 

1.68 

0.00 

0.00 

0.00 

0.00 

0.10 

0.84 

1.88 

1.21 

0.36 

4li 

— 

2.13 

0.06 

1.42 

1.02 

0.00 

0.87 

0.37 

0.14 

0.29 

A 2 

— 

— 

0.01 

0.00 

0.00 

0.00 

0.01 

0.00 

0.00 

0.00 

A 3 

— 

— 

— 

0.82 

0.73 

0.00 

0.13 

0.30 

0.10 

0.15 

A 4 

— 

— 

— 

— 

0.61 

0.00 

0.88 

0.51 

0.89 

0.14 

^5 

— 

— 

— 

— 

— 

0.83 

0.27 

0.16 

0.18 

0.16 

^6 

— 

— 

— 

— 

— 

— 

0.67 

0.00 

0.00 

0.00 

A 7 

— 

— 

— 

— 

— 

— 

— 

0.86 

0.00 

0.00 

As 

— 

— 

— 

— 

— 

— 

— 

— 

1.06 

0.00 


- 

- 

- 

- 

- 

- 

- 

- 

- 

0.18 


Table 29 : Model-dependent systematic uncertainties on the interference fit fractions (%) of the 
resonant contributions for the K-matrix model with m( 74 *^ 7 r^) > 2.1 GeV/c^. The resonances 
are: (^ 0 ) K-matrix S-wave, (^i) / 9 ( 770 ), (^2) w( 782 ), (^43) / 9 ( 1450 ), (^44) / 3 ( 1700 ), (^5) / 2 ( 1270 ), 
{Aq) 74 ° 7 r“ P-wave, (^7) 74 g( 2400 )“, (Mg) 742 ( 2460 )“, (^ 9 ) 743 ( 2760 )“ The diagonal elements 
correspond to the statistical uncertainties on the fit fractions shown in Table 



Ao 

Ai 

A2 

A3 

A4 

As 

As 

A7 

As 

A^ 

A{) 

1.10 

0.00 

0.00 

0.00 

0.00 

0.02 

0.24 

0.25 

0.40 

0.15 

Ai 

— 

0.79 

0.02 

0.41 

0.25 

0.00 

0.26 

0.29 

0.19 

0.05 

A2 

— 

— 

0.02 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

A3 

— 

— 

— 

0.21 

0.08 

0.00 

0.20 

0.09 

0.12 

0.03 

Aa 

— 

— 

— 

— 

0.12 

0.00 

0.14 

0.12 

0.08 

0.02 

^5 

— 

— 

— 

— 

— 

0.58 

0.14 

0.19 

0.08 

0.08 

^6 

— 

— 

— 

— 

— 

— 

0.75 

0.00 

0.00 

0.00 

A7 

— 

— 

— 

— 

— 

— 

— 

0.52 

0.00 

0.00 

As 

— 

— 

— 

— 

— 

— 

— 

— 

0.54 

0.00 

A^ 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0.07 
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D Results of the K-matrix parameters 


The moduli and phases of the K-matrix parameters in Eq. are listed in Table 30 The 


break-down of systematic uncertainties are shown in Tables 31 and 32 


Table 30: The moduli and phases of the K-matrix parameters. The first uncertainty is statistical, 
the second the experimental systematic, and the third the model-dependent systematic. The 
moduli are normalised to that of the p(770) contribution and the phase of p(770) is set to 0°. 


Parameter Modulus Phase (°) 


/lO 

17.0 

± 

3.3 

± 

9.5 

± 

3.7 

347.3 

± 

13.7 

± 

18.7 

± 

3.2 

/ll 

14.9 

± 

17.1 

± 

20.3 

± 

8.0 

160.0 

± 

70.1 

± 

39.6 

± 

26.2 

/l2 

111.3 

± 

23.1 

± 

23.8 

± 

12.8 

226.1 

± 

12.0 

± 

11.2 

± 

4.9 

/l3 

28.7 

± 

14.2 

± 

8.1 

± 

5.3 

186.5 

± 

30.0 

± 

30.4 

± 

8.6 

/l4 

31.0 

± 

12.8 

± 

13.4 

± 

8.8 

10.61 

± 

25.9 

± 

15.2 

± 

2.9 

/3o 

9.5 

± 

1.8 

± 

2.9 

± 

1.1 

20.7 

± 

15.2 

± 

13.5 

± 

10.4 

/3i 

17.2 

± 

6.4 

± 

6.2 

± 

4.8 

19.6 

± 

19.4 

± 

14.4 

± 

3.7 

h 

34.9 

± 

7.6 

± 

14.3 

± 

3.1 

128.3 

± 

12.1 

± 

2.1 

± 

1.9 

A 

53.5 

± 

14.3 

± 

9.2 

± 

4.2 

138.7 

± 

15.5 

± 

7.2 

± 

3.9 

/34 

52.5 

± 

10.2 

± 

22.4 

± 

5.9 

305.0 

± 

10.5 

± 

13.5 

± 

2.2 


Table 31: Systematic uncertainties on the moduli of the K-matrix parameters. The moduli are 
normalised to that of /9(770). 


Source 

fio 

/ll 

fl2 

/l3 

/l4 

do 

di 

d2 

ds 

d4 

PID 

8.1 

5.5 

22.5 

3.2 

13.2 

2.5 

3.6 

10.7 

8.0 

18.0 

Trigger 

0.4 

1.1 

2.5 

1.0 

1.2 

0.2 

0.7 

0.6 

1.1 

1.0 

Reconstruction 

0.1 

0.9 

0.4 

0.5 

0.7 

0.1 

0.4 

0.4 

0.6 

0.5 

Simulation statistic 

0.3 

4.6 

2.8 

2.7 

1.6 

0.1 

1.8 

1.0 

0.6 

0.7 

Background model 

0.6 

0.5 

6.5 

< 0.1 

0.2 

0.1 

0.7 

1.7 

2.4 

2.7 

i4*(2010)- veto 

4.9 

18.9 

1.9 

6.8 

0.6 

1.4 

4.6 

9.2 

3.7 

13.1 

Total (experiment) 

9.5 

20.3 

23.8 

8.1 

13.4 

2.9 

6.2 

14.3 

9.2 

22.4 

Additional resonances 

3.5 

8.0 

12.6 

3.8 

7.7 

1.1 

4.6 

3.0 

3.5 

5.6 

RBW parameters 

1.2 

0.4 

2.3 

3.3 

4.0 

0.2 

1.1 

0.5 

2.2 

1.9 

TTTT res. mass, width 

0.3 

0.7 

0.1 

1.5 

1.7 

0.2 

0.6 

0.4 

0.7 

0.3 

Total (model) 

3.7 

8.0 

12.8 

5.3 

oq 

00 

1.1 

4.8 

3.1 

4.2 

5.9 

Total (all) 

10.2 

21.8 

27.0 

9.6 

16.0 

3.1 

7.8 

14.6 

10.1 

23.2 
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Table 32: Systematic uncertainties on the phases (°) of the K-matrix parameters. The phase of 
p(700) is set to 0° as the reference. 


Source 

fio 

/ii 

/l2 

/l3 

/l4 

/3o 

/3i 

P 2 

/?3 

Pi 

PID 

18.0 

35.0 

9.6 

25.9 

5.8 

10.4 

10.4 

1.7 

1.4 

4.3 

Trigger 

0.5 

0.2 

0.3 

0.2 

< 0.1 

0.5 

0.3 

0.1 

0.5 

0.1 

Reconstruction 

0.3 

2.4 

0.6 

0.5 

0.6 

0.6 

0.4 

0.1 

0.6 

0.2 

Simulation statistic 

4.0 

1.0 

0.2 

2.0 

7.1 

5.2 

4.6 

0.8 

1.9 

1.3 

Background model 

3.1 

8.2 

0.4 

3.9 

0.4 

1.9 

2.5 

0.1 

0.3 

0.4 

D*{2010)- veto 

0.3 

16.4 

5.8 

15.2 

12.1 

6.5 

8.5 

1.0 

6.7 

12.7 

Total (experiment) 

18.7 

39.6 

11.2 

30.4 

15.2 

13.5 

14.4 

2.1 

7.2 

13.5 

Additional resonances 

0.4 

5.0 

4.1 

2.1 

1.0 

3.6 

1.9 

0.2 

3.7 

1.9 

RBW parameters 

3.2 

25.2 

2.6 

8.3 

0.6 

9.2 

3.2 

1.9 

1.0 

1.1 

TTTT res. mass, width 

0.3 

5.0 

0.5 

1.0 

2.6 

3.2 

0.4 

0.2 

0.7 

0.3 

Total (model) 

3.2 

26.2 

4.9 

8.6 

2.9 

10.4 

3.7 

1.9 

3.9 

2.2 

Total (all) 

19.0 

47.5 

12.3 

31.6 

15.5 

17.0 

14.9 

2.9 

8.1 

13.7 
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